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1.  Introduction 


1.1.  Overview 


In  the  absence  of  a  general  analytical  design  tool  for  the  purpose  of  judging 
the  allowable  extent  of  corrosion  wastage  in  oil  tankers,  this  particular  area 
of  ship  structural  analysis  has  been  given  to  experience-based  guidelines  as 
the  only  feasible  treatment.  In  the  rule  books  of  today's  classification 
societies,  the  subject  of  "allowable  wastage"  is  generally  absorbed  into  simple 
equations  that  provide  some  indication  of  a  minimum  strength  standard  for 
newbuild  designs  and  renewals.  While  safe  ships  have  been  built  and  are 
continuing  to  be  built  under  the  guidelines  of  these  rule  requirements,  the 
provisions  involve  a  startlingly  simple  set  of  vEiriables  when  one  considers 
the  complexity  and  diversity  of  the  structures,  the  environments,  and  the 
operation  philosophies  involved  in  today’s  tanker  trade. 

This  report  summarizes  the  work  done  under  the  sponsorship  of 
Ishikawsuima-Harima  Heavy  Industries  (IHI)  and  Mitsubishi  Heavy 
Industries  (MHI)  to  develop  a  rational  approach  to  defining 
corrosion  limits  in  tankers.  The  objective  of  this  project  was  to  make 
advances  in  the  area  of  setting  allowable  limits  for  the  wastage  of 
tanker  structures  based  on  a  procedure  involving  rational  analytical 
techniques  as  an  adjunct  to  the  traditional,  experienced  based 
approaches. 

1.2.  Review  OF  Current  Practice 

The  prediction  of  the  actual  loss  in  the  structural  capacity  of  the  ship, 
structure  due  to  corrosion  can  only  be  dealt  with  on  a  case-by-case  basis.  The 
prediction  must  be  based  on  the  full  facts  of  each  specific  design.  The 
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methods  by  which  these  predictions  are  made  are,  by  nature,  unique  to  each 
design  and  no  attempt  will  be  made  to  review  them.  What  is  reviewed  here 
are  the  criteria  set  out  by  the  classification  societies  which  define  minimum 
requirements  for  hull  strength  and  how  appropriate  forms  of  corrosion 
control  can  result  in  allowances  for  scantlings  below  the  minimum  values. 

A  unified  hull  girder  longitudinal  strength  standard  has  been  established  by 
the  International  Association  of  Classification  Societies  (lACS)  which  all 
ships,  new  and  existing,  must  satisfy.  This  standard,  which  comes  in  the 
form  of  a  simple  formula  for  the  minimvun  midship  section  modulus, 
embodies  the  vast  experience  that  has  been  accmnulated  by  the  members  of 
the  classification  societies  and  has  provided  adequate  safety  for  the  world  s 
fleet  of  ocean-going  vessels.  The  standard  was  most  recently  revised  in  1989 
and  is  as  follows; 

S7  Minimum  Longitudinal  Strength  Standards 

S7.1  The  minimum  midship  section  modulus  at  deck  and  keel  for  ships  90  m  < 
L  <  500  m  and  made  of  hull  structural  steel  is: 

=cL'B(C,+0.7)k  (cm^) 

where  L  =  rule  length  (m) 

B  =  rule  breadth  (m) 

Cb  =  rule  block  coefficient  ( >  0.60  ) 
c  =  Cn  for  new  ships 
c  =  Cj  for  ships  in  service  =  0.9Cn 


=  10.75-| 

o 

o 

1 

for  90m  <  L  ^  300m 

.100  J 

=  10.75 

for 300m  <  L  <  350m 

1 

OJ 

o 

o 

=  10.75- 

for  350m  ^  L  <  500m 

V  100  , 

J 

k  =  material  factor 
=  1.0  for  ordinary  hull  steel 
<1.0  for  higher  tensile  steel 


5 


57.2  Scantlings  of  all  continuous  longitudinal  members  of  hull  girder  based  on 
the  section  modulus  requirement  in  S7.1  are  to  be  maintained  within  0.4L 
amidships. 

However,  in  special  cases,  based  on  consideration  of  type  of  ship,  hull  form  and 
loading  conditions,  the  scantlings  may  be  gradually  reduced  towards  the  end  of 

the  0.4L  part,  bearing  in  mind  the  desire  not  to  inhibit  the  vessel's  loading 
flexibility. 

57.3  In  ships  where  part  of  the  longitudinal  strength  material  in  the  deck  or 
bottom  area  are  forming  boundaries  of  tanks  for  oil  cargoes  or  ballast  water  and 
such  tanks  are  provided  with  an  effective  corrosion  protection  system,  certain 
reductions  in  the  settlings  of  these  boundaries  are  allowed.  These  reductions, 
however,  should  in  ho  case  reduce  the  rhinimum  hull  girder  section  modulus  for  a 
new  ship  by  more  than  5%. 

By  establishing  this  strength  standard  based  on  the  acquired  experience  of 
successful  designs,  a  safety  margin  to  account  for  the  inevitable  wastage  of 
hull  steel  structures  has  been  built  in  to  the  formula.  The  individual 
classification  societies  then  go  on  to  provide  exceptions  to  the  rule  to  account 
for  unusual  design  concepts  and  the  use  of  corrosion  protection  systems  [1] 

1.3.  Problem  Definition 

It  is  clear  that  much  still  needs  to  be  done  to  study  the  problem  that  corrosion 
presents  to  tanker  structures  even  before  solutions  can  be  obtained.  There 
are  many  sources  of  uncertainties  that  are  involved  in  this  particular  aspect 
^sing  of  a  vessel,  and  their  effects,  when  combined,  lead  to  a  very 
complicated  problem.  The  challenge,  therefore,  is  two-fold.  First,  an  overall 
approach  must  be  developed  to  coordinate  the  vast  amount  of  information, 
data,  and  general  theoretical  concepts  involved.  Secondly,  each  component, 
each  module  of  the  procedure  must  be  generated  using  the  most  efficient  and 
accurate  analytical  tools  and  theories  available  given  the  limitations  of 
computational  resources. 

There  is  a  vast  difference  between  the  structural  analysis  of  a  particular 
vessel  under  specific  conditions  and  the  general  treatment  of  an  entire 
tanker  fleet.  Highly  sophisticated  proven  techniques  are  available  to 
accurately  predict  the  strength  of  ships’  structures.  For  example,  non-linear 
finite  element  analyses  exist  to  compute  the  capacity  of  steel  structures  to 
resist  failure  in  a  variety  of  failure  modes,  and  numerical  techniques  are 


6 


7 


available  to  accurately  describe  the  loading  environment  and  load  effects,  but 
these  techniques  are  only  applicable  to  highly  detailed  case-specific  studies, 
and  they  come  at  great  cost  in  computer  time  and  resources. 

In  contrast,  as  mentioned  previously,  the  foimdations  of  the  classification 
society  corrosion  wastage  criteria  and  structural  guidelines  consist  of  very 
broad  general  methods  that  can  only  be  used  as  g^dance.  While  these 
guidelines  provide  a  quick  evaluation  of  a  newbuild’s  performance  or  an 
existing  ship’s  condition,  they  have  no  rational  analytical  basis,  and,  as 
quantified  by  Shama  [2]  a  large  undue  cost  can  be  potentially  developed  as  a 
consequence  of  an  irrationally  designed  structure. 

What  follows  is  a  description  of  the  attempt  made  during  this  one-year 
project  to  bridge  the  gap  between  the  specific  and  general  methods  of 
determining  corrosion  wastage  limits.  This  implies  the  development  of  a 
rational  analytical  tool  that  is  not  too  expensive  to  use,  can  be  used 
interactively  (as  in  the  early  stages  of  design  or  during  a  routine  inspection), 
and  can  be  applied  to  the  general  tanker  fleet.  It  is  with  this  goal  in  mind 
that  the  project  was  undertaken. 

1.4.  Solution  Algorithm 
1.4.1.  Life  Assessment 

The  task  of  defining  corrosion  limits  for  a  complicated  structiiral  system  such 
as  an  oil  tanker  is,  in  fact,  just  one  aspect  of  what  can  be  generally  termed 
the  Life  Assessment  paradigm.  A  particuleu*  vessel  can  age  in  many  ways 
over  its  lifetime,  and  the  purpose  of  a  life  assessment  is  to  develop  some 
global  index  that  describes  the  condition  of  the  aging  vessel  in  terms  of  safety 
or  reliability  or  serviceability,  etc.  It  follows  naturally  that  the  development 
of  a  life  assessment  procedxire  will  provide  a  convenient  framework  from 
which  to  begin  defining  these  corrosion  limits. 

The  main  result  of  a  life  assessment  is  a  description  of  how  the  defined  SI 
(inverse  measure  of  the  probability  of  "failure")  behaves  as  the  vessel  ages, 
i.e.  it  determines  P(t),  where  P  is  the  safety  index  and  is  a  function  of  time. 
Once  P  falls  theoretically  below  a  predetermined  minimum  level,  Pinin>  the 


7 


r 


time  at  which  this  occurs  can  be  noted,  and  limits  can  be  prescribed  based  on 
how  p  was  defined  and  what  caused  it  to  drop  below  the  allowable  level.  It  is 
important  to  note  that  the  initial  limit  (Pmin)  is  determined  firom 
considerations  such  as  economic,  political,  and  social  issues.  It  is  then  the 
variables  that  constitute  the  definition  of  P  in  which  the  engineer  is 
interested  and  to  which  limits  will  be  assigned. 

1.4.2.  Time  Variability  and  Corrosion  Rates 

The  time  dimension  in  this  particular  application  of  life  assessment  methods 
is  constructed  by  the  inverse  of  corrosion  rates  multiplied  by  steel 
thicknesses.  For  this  reason,  accurate  corrosion  rates  are  an  essential  part  of 
this  project.  In  the  first  year  of  the  Structural  Maintenance  for  New  and 
Existing  Ships  Project,  Pollard  focused  on  the  determination  of  corrosion 
rates  in  tanker  internal  structures.  A  large  amount  of  wastage  data  was 
gathered  from  a  wide  range  of  gauging  reports.  Statistical  analyses  were 
performed  to  determine  corrosion  rate  trends  based  on  the  type  of  tank,  the 
type  of  structural  detail,  and  the  relative  location  of  the  detail  within  a  tank. 
It  is  this  type  of  information  that  facilitates  a  more  realistic,  rational  view  at 
monitoring  the  decline  of  a  vessel's  structure  over  the  course  of  its  design 
lifetime. 

1.4.3.  Reliability 

Gauging  corrosion  rates  and  their  effect  on  ships'  structures  is  a  very 
uncertain  proposition,  therefore,  it  is  practically  useless  to  approach  this 
problem  from  a  purely  deterministic  point  of  view.  Any  overall  safety  index 
that  can  be  applied  to  this  problem  is  itself  is  a  random  variable  in  the 
extreme,  and  the  imcertainties  grow  significantly  the  further  into  the  future 
that  P  is  projected.  Much  of  the  uncertainty  is  simply  inherent  in  this  very 
complex  problem  and  can  not  be  reduced.  However,  a  large  portion  of  the 
uncertainty  will  come  firom  modeling  errors  which  reflect  the  limitations  of 
the  available  theories. 

1.5.  Ship  Maintenance  Information  System  (SMIS) 

1.5.1.  Overview 
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The  ever  increasing  availability  of  computer  resources  and  the  growing 
refinement  of  analytical  techniques  make  it  possible  to  take  a  more 
anal3dacal  angle  at  the  problem  of  predicting  that  point  in  a  vessel's  life  when 
the  degree  of  corrosion  wastage  renders  the  structure  unreliable.  Naturally, 
the  design  of  a  computer  application  is  an  integral  part  of  any  attempt  to 
develop  a  solution  to  a  problem  of  such  complexity  as  this  one.  Therefore, 
along  with  the  engineering  considerations  involved  in  this  project,  the 
preliminary  design  of  an  information  system,  the  Ship  Maintenance 
Information  System,  is  interwoven  with  the  theory. 

The  goal  implicit  in  the  development  of  such  a  system  is  to  develop  a  PC  application 
based  on  the  developed  approach  which  has  the  following  qualities: 

•  Efficiency  -  Intelligent  use  of  available  resources. 

•  Flexibility  •  Built  in  capability  for  customizing  the  system. 

•  Reliability  •  Robust  system  with  error  checking  and  input  validation 

procedures. 

•  Maintainability  -  Clear  and  complete  system  documentation  both  of 
the  system  design  and  implementation. 

•  Usability  -  Can  be  applied  by  a  wide  range  of  users. 

•  Accuracy  -  Yields  reasonable  and  useful  restdts. 

As  a  supplement  to  the  theoretical  effort  aimed  at  achieving  the  stated 
objectives,  the  development  of  a  model  program  was  proposed  to  illustrate  the 
point.  This  report,  therefore,  also  documents  the  development  of  the  Ship 
Maintenance  Information  System  (SMIS),  a  PC  based  system  that  was 
modeled  after  the  theoretical  approach  developed  during  this  one  year 
project.  The  SMIS  is  intended  to  be  an  illustration  of  how  such  an  approach 
could  be  implemented. 

1.5.2.  Primary  Programming  Considerations 

The  lifetime  structural  characteristics  of  a  vessel  fleet  constitutes  an 
extremely  complex  physical  situation  which,  to  model,  represents  a 
formidable  and  sometimes  overwhelming  task.  The  amount  of  data  required 
to  represent  even  one  year  of  a  vessel's  life  could  fill  volumes.  In  order  to 
treat  the  many  aspects  of  this  subject,  it  is  necessary  to  break  this  large 
amount  of  data  up  into  small  pieces  that  can  be  handled  one  step  at  a  time  in 
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manageable  portions.  A  relational  database  immediately  suggests  itself  as  a 
means  by  which  to  achieve  this  organization. 

FoxPro  for  Windows  is  a  Microsoft  relational  database  management  system. 
In  addition  to  providing  the  tools  by  which  to  manage  large  amounts  of 
related  data,  FoxPro  also  provides  a  programming  Ismguage  which  allows 
the  development  of  a  sophisticated  user  interface  and  the  precise  control  of 
information  flow.  With  these  powerful  capabilities  provided,  the  entire 
application  could  be  developed  from  within  the  FoxPro  environment. 

However,  while  the  underlying  data  structure  is  easily  constructed  and  the 
management  of  the  data  can  be  framed  in  a  "user  friendly"  interface,  there 
are  a  number  of  aspects  of  the  procedure  that  involve  a  significant  amount  of 
number  crunching,  or  the  repeated  manipulation  of  large  data  sets.  These 
procedures  are  not  suited  to  the  data  management  environment,  but  rather 
to  the  speed  and  simplicity  of  FORTRAN  programming. 

1.5.3.  Design  Limitations 

In  an  attempt  to  design  this  application,  it  is  important  to  realize  the 
limitations  that  are  implicit  in  the  scope  of  this  one  year  project.  Only  the 
first  of  the  two  main  challenges  stated  in  the  overview  was  addressed,  i.e. 
only  the  general  approach  was  modeled.  The  scope  of  the  rigorous  technical 
aspects  was  reduced  to  ensure  that  the  design  itself  was  completed.  In  view 
of  this,  the  following  general  simplifications  were  made: 

•  It  was  not  possible  to  address  all  of  the  failure  modes  that  are  the  result  of 
corrosion  in  hull  structures.  The  strength  (capacity)  analyses  were 
focused  on  failure  due  to  buckling  instability  of  the  ships’  structural 
components.  Failure  due  to  corrosion  fatigue  and  cracking  were  not  dealt 
with  directly. 

•  The  treatment  of  corrosion  rates  was  limited  to  general  uniform  wastage. 
Pitting  and  grooving  types  of  corrosion  were  not  treated 

•  Simplified  Reliability  Methods  were  used  to  limit  the  complexity  of  the 
System  Reliability  problem  to  a  manageable  level. 
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1.6.4.  Context  Diagram 

Shown  in  Figure  1.1  is  the  context  layer  diagram  for  the  theoretical  SMIS 
application.  A  context  layer  data  flow  diagram  represents  the  general 
interface  between  the  application  and  the  external  sources  and  sinks  of 
information.  This  particular  diagram  shows  that  a  combination  of  vessel 
specific  data  and  fleet  wide  data  are  input  into  the  system  where  they  are 
used  to  describe  the  availability  (a  general  description  of  reliability)  of  the 
vessel  projected  over  time.  The  system  then  generates  a  report  of  corrosion 
limits  based  on  the  results  of  the  life  assessment.  The  components  of  this 
system  will  be  developed  over  the  next  three  chapters. 


Figure  1.1:  SMIS  context  layer  diagram 
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2.  Life  Assessment 


2.1.  Availabiuty 

A  life  assessment  procedure  provides  a  convenient  framework  from  which  to 
prescribe  limiting  conditions  on  any  one  of  the  many  factors  that  affect  the 
reliability  of  a  vessel.  Nippon  Kaiji-Kyokai  developed  a  model  of  a  life 
assement  procedure  for  ships  and  offshore  structures  that  could  theoretically 
provide  a  comprehensive  indication  of  the  condition  of  a  particular  vessel  at 
any  one  time  dining  its  operational  lifetime.  In  this  life  assessment 
approach,  the  reliability  is  defined  in  terms  of  the  availability  of  the  vessel, 
a  requirement  set  by  the  owners/operators  that  describes  the  percentage  of 
time  that  the  vessel  must  be  in  service. 

During  a  ship's  lifetime,  it  spends  a  certain  amount  of  time  being  inspected 
or  repaired.  These  ’outages"  can  be  attributed  to  three  major  categories  of 
events: 

1.  Planned  Inspection  and  Maintenance  Routines  (IMR)  either  required  by 
law  or  set  by  the  owners  themselves  (whichever  is  the  more  conservative 
practice). 

2.  The  repair  of  structural  failures  that  are  due  to  a  weakness  in  the  ship's 
structure.  These  outages  become  more  frequent  as  the  ship  ages. 

3.  The  repair  of  structural  failures  following  accidents  that  are  caused  by 
unforeseen  extreme  -loading  conditions  and/or  human  and  organizational 
error  (HOE). 

A  numerical  quantity  called  the  unavailability  can  be  defined  as  that  fraction 
of  time  that  the  vessel  is  out  of  service  (years-per-year)  due  to  each  of  the 
above  three  categories.  Respectively,  these  components  of  the  total 
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unavailability,  U,  can  be  designated  as  t/p^,  Ugf,  Uqt-  The  availability.  Ay,  is 
expressed  as: 


Av  ss  1  - 17  ss  1 .  (Upl  +  Ugp  +  Uqj-) 

If  a  design  Av  is  given,  and  provided  that  the  components  of  unavailability 
can  be  accurately  calculated  or  predicted  over  the  life  of  the  vessel, 
judgments  can  be  made  concerning  the  acceptable  or  allowable  deterioration 
of  the  vessel’s  structural  strength.  The  figure  below  schematically  shows  this 
process  in  terms  of  the  above  quantities. 


Av  (yrs/yr) 


Figure  2.1  Availability  as  a  Function  of  Time 


In  order  to  chart  the  values  of  unavailability  over  time,  a  combination  of 
detailed  structural  analysis,  experience,  and  a  wealth  of  data  are  needed. 
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2.2.  Support  Database 


A  database  structure  is  needed  to  support  the  types  of  analyses  involved  in 
the  assessment.  The  following  three  major  database  components  serve  as  a 
starting  point  for  the  design  of  the  required  database  struct\ire: 

•  A  preliminary  survey  database  that  would  contain,  among  other  things, 
information  concerning  the  vessels  particulars,  its  cargo,  its  route,  its 
corrosion  protection  system,  its  inspection  and  maintenance  routine,  its 
intended  service  life,  and  its  prescribed  availability,  (design  Av,  Upj) 

•  A  database  of  records  2md  statistics  of  unforeseen  accidents,  instances  of 
human  error  resulting  in  accidents,  etc.  {Uqj) 

•  A  database  containing  referential  data  such  as  gauging  reports,  crack 
inspections,  the  location  and  nature  of  structural  failures,  the  time  it  took 
to  repair  them,  etc.  {Ugf) 

The  nature  of  the  analytical  tool  being  proposed  requires  that  a  database 
management  system  be  designed  to  maintain  the  data  and  control  the  flow  of 
information.  Without  such  a  system,  the  tool  would  be  difiicult  to  employ, 
and  then  only  by  a  small  range  users.  Shown  in  fig.  2.2  is  a  data  flow 
diagram  (DFD)  depicting  the  role  of  the  database  management  system  within 
the  context  of  this  project. 
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Figure  2.2i  First  layer  DFD 


As  indicated  in  fig.  2.2,  the  SMIS  database  management  system  must  be 
designed  to  accept  input  from  a  range  of  users,  allow  an  engineer  to  control 
an  analytical  session,  maintain  and  manage  the  data,  act  as  a  driver  for  the 
analytical  routines,  and  produce  reports  to  ease  the  interpretation  of  the 
results. 
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Shown  in  fig.  2.3  is  a  more  detailed  view  of  the  components  of  the  database 
management  system  that  are  required  to  achieve  its  purpose.  The  exact 
structure  of  the  support  database,  including  the  format  of  the  data  and  how 
it  will  be  used  in  the  analysis  will  be  discussed  in  the  following  chapters. 

2.3.  Un-Availabiuty 

Figure  2.4  shows  the  relationship  between  the  support  database,  the  general 
analysis  modules,  and  the  three  components  of  unavailability. 


rjTn 
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2.3.1.  UpL:  Planned  Outages 


The  unavailability  due  to  planned  outages  and  the  vessel's  IMR  can  be 
derived  from  the  information  that  is  contained  in  the  Preliminary  Survey 
database.  This  quantity  which  may  vary  with  time  (the  owner  might  think  it 
necessary  to  decrease  the  amo\mt  of  time  between  inspections  as  the  ship 
gets  older)  must  be  supplied  to  the  database. 

2.3.2.  Uqt*  Human  Error  and  Other  Causes 

The  unavailability  due  to  accidents  and  human  error  can  be  evaluated  based 
on  past  experience.  This  is  the  reason  for  the  database  containing  records 
and  statistics  of  such  events. 

2.3.3.  Ugp  :  Unavailability  due  to  Structural  Failure 

The  majority  of  the  analytical  effort  undertaken  in  this  project  surrovmded 
the  calculation  of  the  unavailability  due  to  structural  failure.  This  effort 
involves:  collecting  and  categorizing  the  incidents  of  failure,  providing  a 
statistical  interpretation  of  the  corrosion  wastage  data,  developing  a 
statistical  model  of  the  prescribed  loading  condition,  developing  a  best 
estimate  of  the  structimal  capacity,  and  finally,  through  reliability  methods, 
obtaining  the  annual  probability  of  failure  for  each  mode  of  failure  and  for 
each  year  of  the  service  life. 

Ugf,  as  defined  by  NK,  is  as  follows: 

u  u)-i  _ 

S  MTTR,  +  MTBF,  (0 

where  there  are  7  failure  modes,  and  MTTR  and  MTBF  are  defined  by  the 
following: 

MTTRi:  Mean  time  to  repair  failure  i.  (obtained  from  the  structural 
failure  incident  database) 

MTBFi :  Mean  time  between  failures  in  the  i^  mode. 
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The  mean  time  between  failure  for  a  parti cvilar  mode  (in  years)  is  simply  the 
inverse  of  the  anual  probability  of  failure  for  that  mode.  The  calculation  of 
the  probability  of  failure  in  a  particular  mode  is  a  subject  of  reliability 
analysis  that  is  well  known  throughout  the  industry.  For  a  given  loading 
condition  (Demand)  and  a  predicted  structural  strength  (Capacity)  there  are 
several  levels  of  complexity  that  may  be  employed  to  obtain  the  probability  of 
failure.  There  is  a  great  deal  of  analysis,  judgment,  and  experience  that  is 
required  before  meaningful  results  can  be  realized. 

2.4.  Failure  Modes 

For  the  purpose  of  estimating  Usp,  it  is  necessary  to  collect  structural  failure 
incidences  into  general  categories  from  which  information  can  be  drawn  that 
will  be  applicable  to  any  vessel  in  the  fleet.  In  reality,  no  two  structural 
failure  incidences  are  exactly  alike.  However,  these  incidences  can  be 
classed,  and  it  is  these  classes  or  modes  of  failure  upon  which  the  analytical 
tool  will  operate.  According  to  Daidola  et.  al.  [3]  in  terms  of  the  longitudinal 
stregth  of  a  hull  girder,  there  are  five  general  categories  of  failure: 

•  Yield  failure  due  to  bending  of  the  ship  considered  as  a  beam 

•  Compression  instability  buckling 

•  Brittle  fracture 

•  Fatigue  fracture 

•  Ultimate  plastic  collapse 


These  five  general  modes  can  be  further  separated  into  categories  based 
primarily  on  the  type  of  structural  sub-elements  that  sire  affected.  Only 
compression  instability  buckling  and  laltimate  collapse  are  treated  in  this 
study;  although,  given  the  appropriate  support  data  and  analytical 
techniques,  the  overall  procedure  could  be  extended  to  treat  the  other  three 
general  categories  of  failure. 
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It  is  necessary  to  calculate  the  Mean  Time  Between  Failure  (MTBF)  and  the 
Mean  Time  To  Repair  (MTTR)  for  each  failure  mode  that  can  be  identified  as 
being  a  likely  to  occur  during  the  design  lifetime  of  the  vessel.  The 
development  of  MTBF  for  a  particular  mode  is  an  anal3rtical  matter  that  will 
be  discussed  in  the  next  chapter.  MTTR,  on  the  other  hand,  must  be 
obtained  from  fleet- wide  data  and  experience  in  a  manner  similar  to  that  for 
obtaining  corrosion  rate  information.  Therefore,  under  the  heading  of 
Referenda!  Data,  the  support  database  must  contain  information  in  a  form 
that  will  yield  appropriately  categorized  repair  information 

This  categorization  process  requires  a  great  deal  of  shipyard  experience  and 
data  and  could  potentially  be  carried  out  to  a  high  level  of  detail. 
Unfortunately,  this  type  of  data  is  generally  held  as  confidential  and  was  not 
available  during  this  study.  However,  some  additional  general  categorization 
can  be  made  which  can  serve  as  an  illustration  and  a  starting  point  for 
further  work  on  this  topic. 

In  the  case  of  compressive  instability  buckling,  repair  information  can  be 
seperated  into  the  following  five  general  categories: 

Class  I  -  failure  leading  to  the  replacement  of  longitudinal  stiffeners 
(tripping,  stiffener  induced  buckling,  plate  induced  buckling) 

Class  II  -  failure  leading  to  the  replacement  of  internal  plating 
between  stiffeners  (buckling  of  plating  between  stiffeners) 

Class  III  -  failure  leading  to  the  replacement  of  shell  (external)  plating 
between  stiffeners  (buckling  of  plating  between  stiffeners) 

Class  IV  -  failure  leading  to  the  replacement  of  an  internal  stiffened 
panel  (overall  grillage  buckling) 

Class  V  -  failure  leading  to  the  replacement  of  an  external  stiffened 
panel  (overall  grillage  buckling) 

While  these  five  classes  cover  nearly  all  of  the  types  of  compressive 
instability  buckling  failures,  some  additional  information  must  be  supplied 
in  order  to  get  accurate  information  regarding  how  much  time  a  particular 
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vessel  will  be  out  of  service  because  of  them.  For  example,  there  is  a  high 
degree  of  correlation  between  failure  modes  and  in  the  case  of  an  entire 
stiffened  panel  being  replaced,  including  the  time  it  takes  to  repair  each 
stiffener  on  the  buckled  panel  would  result  in  an  overestimation  of  the  repair 
time.  These  difficulties  can  be  dealt  with  but  only  if  the  required  information 
is  provided. 
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3.  Structubal  Failure,  Usf 


3.1.  Overview 

The  determination  of  the  unavailability  due  to  year-to-year  type  structural 
failures,  Usp,  comprises  the  major  analytical  effort  of  this  Life  Assessment 
routine.  A  large  amount  of  data  analysis  as  well  as  theoretical  concepts  are 
required  to  model  a  particular  vessel’s  service  lifetime. 

As  mentioned  in  the  previous  chapter,  Ugp  has  been  defined  by  NK  as  a 
function  of  the  mean  time  between  failure  incidences  and  the  mean  time  that 
the  vessel  is  unavailable  while  the  failure  is  being  repaired. 

Since  specific  types  of  failure  tend  to  occur  more  often  as  a  vessel  ages,  UgF 
is  a  function  of  time  and  the  particular  failure  modes  that  are  associated  with 
the  vessel.  MTTRn  (Mean  Time  To  Repair  failure  mode  'n')  is  a  quantity  that 
is  obtained  for  each  failure  mode  through  the  analysis  of  Inspection  and 
Maintenance  Routine  (IMR)  data  collected  and  stored  in  the  support 
database  and  will  be  assumed  to  be  constant  over  the  life  of  the  vessel  being 
examined.  MTBFn  (Mean  Time  Between  Failure  mode  'n')  is  cast  in  terms  of 
years,  and  is  defined  as  the  inverse  of  the  annual  propability  of  failure  for 
the  n^^  failure  mode.  These  probabilities  will  increase  in  time  due  to 
wastage  of  the  internal  structure.  Since  the  above  summation  is  over  all 
possible  failure  modes,  there  will  be  a  marked  increase  in  Ugp  over  the 
lifetime  of  the  ship. 

3.2.  General  Procedure 

The  task  of  developing  an  estimate  of  the  annual  probability  of  failure  for 
any  given  failure  mode  can  be  divided  into  a  number  of  modules  or 
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subroutines  each  of  which  constitute  a  major  component  of  the  calculation. 
These  modules  are  listed  below: 

•  Vessel  Definition  Module 

•  Failure  Definition  Module 

•  Capacity  Module 

•  Demand  Module 

•  Reliability  Module 

•  Corrosion  Module 


The  general  procedure  involves  defining  a  section  of  a  particular  vessel’s  hull 
and  the  failure  modes  associated  with  it.  For  example,  many  incidences  of 
buckling  occur  around  the  midship  section  where  the  primary  bending 
moment  is  generally  at  its  peak.  These  incidences  might  range  fi'om  very 
localized  buckling  of  plating  between  stiffeners  to  overall  collapse  of  the 
primary  structure. 

Next,  the  loading  effects  are  determined  based  on  a  particular  vessels 
geometry  and  loading  environment.  The  capacity  of  the  structure  and  its 
elements  are  then  calculated  and  compared  with  the  demands  of  the  seaway 
loads.  This  involves  the  use  of  reliability  methods  that  treat  both  individual 
structural  elements  and  systems  of  elements.  Combining  knowledge  of  the 
resulting  probability  of  failiu*e  and  knowledge  of  the  consequences  (repair 
time)  of  failure  for  each  mode  results  in  a  calculation  of  llgp  for  one  given 
time  step. 

Using  the  corrosion  data  contained  in  the  referencial  database,  corrosion 
rates  can  be  calculated  and  applied  to  each  element  of  the  defined  section. 
The  designated  time  step  defines  the  extent  of  the  wastage  of  these  elements 
and  their  capacities  are  then  recaculated.  Applying  the  same  loads  as  before, 
the  procedure  for  determining  Ugp  is  repeated.  The  next  time  step  is  made, 
the  section  is  corroded  further,  and  the  entire  process  is  repeated  until  Ugp 
is  defined  over  the  entire  Design  Life  of  the  vessel. 
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Figure  3.1  contains  a  data  flow  diagram  for  the  calculation  of  Ugp. 
Components  of  each  module  in  the  procedure  and  their  relationship  of  the 
support  database  components  are  depicted. 

With  the  general  procedure  outlined  above,  what  follows  is  a  description  of 
each  module  and  how  it  fits  in  to  the  calculation  of  UgpCt). 

3.3.  Vessel  Description  Module 

Purpose:  To  provide  all  of  the  vessel  specific  information  that  will  be  needed  as 
input  to  the  analyses  that  follow  in  subsequent  modules. 

In  order  to  complete  the  analyses  outlined  in  the  preceding  section  for  a 
specific  vessel,  the  preparation  of  a  large  amount  of  preliminary  data  is 
required.  Specifically,  extensive  information  must  be  provided  on  the 
physical  structure  of  the  vessel  as  well  as  its  intended  operational 
performance  or  mission  profile. 

The  Physical  Vessel 

There  are  two  main  aspects  involved  in  the  physical  description  of  a  vessel. 
One  aspect  involves  the  description  of  the  hull  geometry  and  weight 
distribution  for  the  pxirpose  of  calculating  Stillwater  and  vertical  wave 
bending  moments,  and  the  other  involves  a  description  of  the  internal 
longitudinal  structural  components  for  the  purpose  of  calculating  the 
capacity  of  the  hull  to  resist  these  moments. 

With  the  obvious  exception  of  the  outside  hull  form,  an  oil  tanker  is  generally 
made  up  of  rectangular  cells.  It  is  divided  internally  by  decks,  transverse 
bulkheads,  and  longitudinal  bulkheads,  which  constitute  planar  divisions 
parallel  to  the  base  plane,  section  plane,  and  centerline  plane,  respectively. 
Therefore,  a  logical  point  to  begin  the  vessel  description  is  with  the 
designation  of  these  imgor  internal  divisions.  The  configuration  of  the  cell 
spaces  is  naturally  a  complicated  one,  and  therefore,  simply  stating  the 
number  of  each  type  of  division  will  generally  not  lead  to  an  accurate 
description  of  the  internal  spaces.  In  order  to  make  this  description  possible, 
it  is  necessary  to  assume  that  a  hierarchy  exists,  i.e.,  one  type  of  structiu'e 
represents  the  primary  division,  another  type  constitutes  the  secondary,  etc. 
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Shown  in  Fig.  3.2  is  one  possible  model  for  this  hierarchy  depicted  as  ”one-to* 
many"  relationships.  In  words  it  states  that  for  every  vessel  there  are  many 
sections  separated  by  transverse  bulkheads,  and  for  every  section  there  are  a 
number  of  deck  levels  separated  by  decks  and  inner  bottoms,  and,  finally,  for 
every  deck  level  there  are  a  number  of  transverse  compartments  separated 
by  longitudinal  bulkheads.  There  will  natiirally  be  configurations  that  can 
not  be  described  by  this  model,  however,  it  is  simple  enough  to  facilitate  a 
quick  and  fairly  realistic  description  of  the  internal  arrangement  of  a  vessel. 


Figure  3.2;  Description  of  the  Internal  Arrangement 


3.3.1.  Description  of  Vessel  for  Load  Calculation 

The  loading  conditions  that  will  be  experienced  by  a  vessel  during  its  lifetime 
are  based  on  the  superposition  of  the  Stillwater  loads  and  the  loads  that  are  a 
result  of  the  vessel's  response  to  its  wave  environment.  For  the  wave  loads  in 
this  study,  only  the  vertical  wave  bending  moment  will  be  examined  as  this 
is  the  primary  component  of  the  axial  stresses  that  cause  buckling  in  a 
vessel's  longitudinal  members.  Other  loads  such  as  transverse  moments  and 
slamming  will  not  be  treated,  although  there  is  room  for  such  analysis  within 
the  overall  procedure. 
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As  mentioned,  the  hull  is  subdivided  longitudinally  into  stations.  Typically, 
there  are  around  20  stations  defined  between  perpendiculars  and  the  result 
is  a  longitudinally  "discretized"  vessel  as  shown  in  fig.  3.3: 
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For  the  calculation  of  both  the  Stillwater  and  vertical  wave  bending  moment, 
the  weight  and  hull  form  are  needed  at  each  station. 
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The  computation  of  the  Stillwater  loads  is  simply  an  element  of  the  basic 
hydrostatic  calculations  that  are  performed  in  every  design  process.  It 
requires  a  knowledge  of  the  longitudinal  weight  distribution  in  a  variety  of 
operating  conditions  (most  importantly;  full  load  and  ballast  conditions)  and 
also  the  outside  form  of  the  "wetted "  hull  in  each  operating  condition  for  the 
purpose  of  calculating  the  bouyancy  distribution. 

3.3.2.  Description  of  Vessel  for  Capacity  Calculation 

A  mathematical  idealization  of  a  ship's  structure  can  be  acheived  in  many 
ways  and  to  many  degrees  of  complexity.  For  the  purpose  of  this  study, 
analyses  are  performed  on  a  single  transverse  cross-section  of  the  vessel  hull 
at  a  time.  This  two-dimensional  structural  model  is  extended  to  three 
dimensions  by  assuming  a  parallel  prismatic  form  between  a  specified 
transverse  web-frame  spacing. 

An  idealized  transverse  section  can  be  subdivided  into  elements  and  groups 
of  elements  whose  structural  response  can  be  estimated  using  established 
theories  and  structural  analysis  techniques.  Combining  system  reliability 
methods  with  these  element  response  analyses  will  lead  to  a  fairly 
comprehensive  treatment  of  a  parallel  section  of  the  hull  from  individual 
panel  buckling  up  to  the  collapse  of  the  primary  structure.  What  follows  is  a 
description  of  a  method  that  can  be  implemented  in  such  an  application. 

Section  Idealization 

Keeping  in  mind  the  assumptions  and  limitations  of  the  structural  analysis 
techniques  to  be  used,  subdiving  one  of  a  vessel's  transverse  sections  requires 
some  judgement  in  order  to  ensure  that  the  the  structural  response 
(buckling)  of  the  resulting  elements  are  accurately  described  by  the  theories. 
The  basic  building  blocks  of  a  longitudinally  fi*amed  parallel  section  of  a 
vessel  consist  of  a  panel  of  shell  plating  along  with  an  attached  longitudinal 
stiffening  structure.  The  term  "element"  used  in  the  context  of  this  study 
applies  to  these  bmlding  blocks  and  examples  are  shown  in  the  figures  below. 
Fig.  3.5  &  3.6  shows  a  cut  out  panel  section  that  could,  for  example,  have 
been  taken  fi-om  the  side  shell.  The  shaded  portion  constitutes  an  "element " 
as  described  above  and,  in  this  particular  figure,  is  representative  of  an 
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6leiiient  thst  h&s  an  L  shapad  longitudinal.  Figure  3.6  shows  the  four  cross 
sectional  configurations  that  will  be  considered  in  this  study. 


Figure  3.5;  Schematic  representafion  nf  an  examnle  element** 


1 .  Plate  with  no  stiffener: 

2.  Plate  with  flat  bar  stiffener: 

1 

3.  Plate  with  T-bar  stiffener: 

T 

4.  Plate  with  'L'-bar  stiffener: 
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Fi£ure  3.6:  Possible  configurations  for  element  cross  section. 


The  structural  analysis  routines  that  will  be  described  require  a  fair  amount 
of  information  about  each  element.  An  individual  element's  dimensions, 
location,  orientation,  and  boundary  conditions  are  all  necessary  ingredients 
and  must  be  accurately  described. 

Dimensions:  An  element's  dimensions  consist  of  the  cross  sectional  shape 
and  plate  thicknesses,  and  the  (longitudinal)  web  frame  spacing.  Shown 
below  is  an  example  of  the  conventions  used  in  this  study: 


Flange  Breadth 


Flange  Thicknessz 


Plate 


Web  Thickness 
Node 


Thickness 


.V 


Web 

Depth 


Plate  Width 


Orientation 
Angle,  Theta 


Figure  3.7:  Element  dimension  conventions 


Location:  An  element's  location  is  defined  as  the  location  of  the  element 
node  with  respect  to  a  coordinate  system  whose  origin  is  defined  as  the 
intersection  of  the  centerline  and  baseline  of  the  section.  An  element's  node 
is  taken  as  the  center  of  the  element's  plate  component  (see  figure  above). 
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Figure  3.8;  Element  location. 


Boundary  Conditions:  An  element's  structural  response  is  strongly 
influenced  by  the  conditions  that  exist  at  its  boundaries.  Through  careful 
definition  of  these  element  boundary  conditions,  it  is  possible  to  model 
element-to-element  interactions  as  well  as  the  presence  of  lateral  loads 
arising  from  hydrostatic  and  internal  cargo  pressure. 

3.3.3.  Mission  Profile 

The  mission  profile  of  a  vessel  outlines  various  information  regarding  the 
vessel’s  operation  requirements,  limitations,  and  expectations.  For  the 
purpose  of  this  study,  the  following  information  is  required: 

•  Design  Lifetime  (yeeirs) 

•  %  of  time  in  Ballast  voyages 

%  of  time  in  Full  Load  voyages 

•  Ballast  Route  (Marsden  Squares  and  associated  time  factor) 

Full  Load  Route  (Marsden  Squares  and  associated  time  factor) 
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•  Operating  Policy:  Speed  vs.  Significant  Wave  Height  for  both 
Ballast  and  Full  Load  Conditions 

•  Inspection  and  Maintenance  Routine 


Shown  below  is  a  schematic  example  of  how  the  planned  unavailability  can 
be  prescribed.  It  reflects  the  possibility  that  planned  inspection  and 
maintenance  might  step  up  over  time. 


Upl  (yrs/yr) 


- - - [-> 

5  10  15  (yrs) 

Fi2ure  3.9:  Planned  unavailability 
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3.4.  Demand  Module 


Purpose:  To  develop  a  probabilistic  model  of  the  extreme  vertical  bending  moment 
for  a  specific  vessel. 

3.4.1.  Overview 

The  "demand"  that  is  imposed  on  a  tanker  vessel  is  made  up  of  many 
different  loading  effects.  Slamming  loads,  Stillwater  bending  moments,  wave 
loads,  and  inertial  forces  all  contribute  to  the  tj^ical  global  loading 
conditions  experienced  by  a  vessel.  In  view  of  longitudinal  strength,  which  is 
the  focus  of  this  study,  only  vertical  bending  moment  will  be  considered  since 
it  contstitutes  nearly  all  of  the  demand  that  is  placed  on  the  longitudinal 
structural  components. 

The  two  principal  components  of  this  vertical  bending  moment  are  the 
StillWater  Bending  Moment  (SWBM)  and  the  Vertical  Wave  Bending 
Moment  (VWBM).  In  deterministic  terms,  the  Total  Vertical  Bending 
Moment  (TVBM)  can  be  expressed  as; 

TVBM  =  SWBM  ±  VWMB 

The  Convention  used  here  is  that  a  negative  value  indicates  a  "sagging" 
moment,  while  a  positive  value  indicates  a  "hogging"  moment.  In  addition,  a 
particular  vessel  is  assumed  to  experience  the  VWBM  symmetrically  in  the 
hogging  and  sagging  mode  (hence,  the  ±  in  the  equation). 

In  reality,  there  is  a  great  deal  of  uncertainty  associated  with  the  above 
equation.  Among  the  many  factors  contributing  to  this  imcertainty,  those 
associated  with  the  inherent  randomness  of  the  ocean  environment  are 
dominant  along  with  the  modeling  errors  that  are  introduced  as  a  result  of 
the  assumption  that  a  ship  responds  linearly  to  its  environment.  In  any  case, 
the  only  rational  approach  to  modelling  the  total  vertical  bending  moment  is 
to  represent  all  of  the  factors  contributing  to  TVBM  in  a  probabilistic  sense 
rather  than  an  exact  mathematic  (deterministic)  sense. 

More  specifically,  both  SWBM  and  VWBM  are  random  variables  and 
therefore,  so  is  TVBM.  The  purpose  of  this  module  therefore  is  to  develop  an 
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expression  for  the  probability  distribution  of  TVBM  for  a  specific  vessel, 
given  its  route  and  response  characteristics.  This  probabilistic 
representation  of  TVBM  (demand)  will  then  be  compared  to  a  similar 
representation  of  the  Capacity  of  the  structure  to  determine  the  failure 
probability. 

Due  to  conservative  design  philosophies,  instability  of  tanker  structures  in 
the  buckling  mode  is  generally  brought  about  only  by  extreme  environmental 
(wave)  conditions.  While  the  Stillwater  loads  can  be  controlled  and 
minimized  to  a  certain  extent,  the  extreme  sea  conditions  make  the  vertical 
wave  bending  moment  the  dominant  load  effect  and  therefore  drives  the 
analysis  of  the  longitudinal  structure.  Much  work  has  been  done  to  develop 
probabilistic  models  of  extreme  sea  conditions  and  their  effect  on  a  vessel, 
and  the  approach  taken  in  this  study  is  based  on  that  developed  by  Mansour 
[4] 


3.4.2.  Environment 

The  first  step  in  this  process  involves  determining  what  sea  conditions  a 
tanker  is  likely  to  face  based  on  available  sea  data.  A  vessel’s  trade  route  can 
be  separated  into  areas  over  which  the  sea  conditions,  typically  charecterized 
by  significant  wave  height,  are  relatively  constant.  There  have  been  a 
number  of  attempts  to  gather  comprehensive  ocean  data,  but  there  has  yet  to 
be  produced  an  adequate  set  of  consistent,  complete  measxirements  from 
which  directional  wave  energy  spectra  can  be  derived.  The  most 
comprehensive  collection  of  measurements  to  date  is  that  compiled  by 
Hogben  and  Lumb  dunng  a  period  of  seven  years  fi*om  1953  to  1961.  Data 
involving  wave  height  and  periods  were  collected  for  areas  that  were  grouped 
into  Mgirsden  square  zones  (shown  below).  An  example  of  their  data  is 
presented  below  for  the  case  of  the  Norwegian  Sea  area  (Marsden  square  #1). 
In  effect,  the  table  represents  a  scatter  diagram  (observed  percentage 
frequency  of  occurance)  of  a  combination  of  wave  height  and  period. 

Table  1.:  Scatter  Diagram  for  Northern  North  Atlantic  Trade  Zone  (Marsden  Squares 
1,2,6,7,&8) 

Wave  Period 
(seconds) 
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With  a  vessel's  mission  profile  outlined  in  the  Vessel  Description  Module,  it 
is  then  possible  to  calculate  the  total  relative  frequency  of  occurance  for  each 
combination  of  significant  wave  height  and  zero  up-crossing  period  based  on 
the  designated  Marsden  Squares  and  the  relative  time  spent  in  each  one  in 
either  the  ballast  and  full  load  conditions.  That  is: 

)  =  £/), («„r.)/(i) 

I 
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where:  pj  =  observed  frequency  of  occurance  of  the  combination  of  Hj 
and  in  Marsden  Sqaure,  /. 

f(i)  =  the  time  factor  (percentage  of  time)  that  the  vessel  spends  in 
Square  i 

The  summation  is  taken  over  all  Marsden  Squares  along  the 
vessel's  route 

A  separate  p(Hs,T2)  matrix  is  formed  for  each  of  the  two  load 
conditions. 

In  order  to  calculate  the  response  of  a  particular  vessel,  each  sea  state  on  the 
vessel  s  trade  route  must  first  be  described  in  terms  of  a  characteristic  wave 
energy  spectrum.  While  the  set  of  wave  records  presented  by  Hogben  and 
Lumb  does  not  provide  enough  information  to  develop  fully  directional  sea 
spectra  for  each  Marsden  square,  there  are  other  idealized  point  spectra  that 
can  be  calculated  from  the  data  and  that  can  provide  valuable  input  to  the 
ship  response  "black  box"  that  will  be  discussed  in  the  next  section. 

Of  the  various  point  spectra  that  are  well  known  to  the  field,  the 
Bretschneider  Spectrum  is  chosen  for  use  in  this  study  since  its  two 
parmeters  (wave  height  and  period)  allow  a  more  accurate  description  of  a 
seaway  than  a  one  parameter  spectrum  (Pierson-Moskowitz)  while  the  sea 
data  available  is  insufficient  for  the  development  of,  say,  the  Ochi  6- 
parameter  spectrum. 

The  Bretschneider  Spectrum  has  the  form  : 

^  0)  Leo  J 

where  the  parameters  A  and  B  are  in  fact  dependent  on  the  parameters  of 
wave  height  and  period.  A  and  B  have  several  forms  depending  on  what 
characteristic  values  for  height  and  period  are  used.  For  example  ,  if  Hogben 
and  Lumb's  data  is  presented  in  terms  of  Significant  wave  height  (Hg)  and 
Zero  Up-crossing  period  (Tj),  then  it  is  convenient  to  express  A  and  B  in 
terms  of  Hg  and  T2  as  follows: 
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A 


243-^\  and  B 

t: 


965 


t: 


The  implication  of  using  an  idealized  point  spectrum  is  that  the  wave  field  in 
the  open  ocean  consists  of  two  dimensional  long  crested  waves.  This  is 
obviously  a  misrepresentation  of  the  real  situation  especially  in  storm 
(extreme)  conditions  and  use  of  this  model  could  potentially  lead  to  a 
significant  overestimation  of  the  environment  to  which  a  vessel  is  subjected. 
The  sea  is  generally  "softened"  by  its  directionality,  and  this  effect  may  be 
partially  accounted  for  by  the  use  of  a  spreading  function  in  conjimction  with 
a  point  spectrum.  At  the  15th  International  Towing  Tabk  Conference  (ITTC 
1978),  it  was  proposed  that  the  Bretschneider  spectrum  be  combined  with  a 


2 

spreading  function  of  the  form:  —  cos'p  (where  to  model 

average  conditions.  The  final  form  of  the  characteristic  wave  energy 
spectrum  for  a  particular  Marsden  Square  is: 


2 

S,  (co,  |i)  = — S,  (cd)  cos^  p 
'  n  ' 

where  S((o)  is  as  defined  previously. 

3.4.3.  Environmental  Effects  (Load) 

With  a  spectral  representation  of  any  given  seaway  established  according  to 
the  preceeding  section,  it  is  possible  to  calculate  a  variety  of  ship  response 
sepctra  for  a  specific  vessel  provided  that  a  Transfer  Fimction  or  Response 
Amplitude  Operator  (R.A.O.)  can  be  developed  for  the  responses  of  the 
particular  vessel.  For  this  study  of  course,  it  is  the  vertical  wave  bending 
moment  at  a  transverse  section  that  is  of  interest. 

The  calculation  of  the  vertical  wave  bending  moment  response  of  a  vessel  at  a 
particular  section  involves  (first)  the  solution  of  the  equations  of  motion  for  a 
ship  in  regular  seas,  (second)  the  evaluation  of  incremental  vertical  forces 
(excluding  Stillwater  buoyant  forces)  based  on  these  motions,  and  (third)  the 
integration  of  these  forces  over  the  length  of  the  vessel. 
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The  vertical  wave  bending  moment  at  a  particular  section  is  equal  to  the 
difference  between  the  inertial  force  and  the  sum  of  the  external  forces: 
exciting  force  (E),  restoring  force,  (R)  and  body  motion  force  (D)  [7]. 

VWBMix)  =  /j  -  (£j  +  /?j  +  D, ) 

Employing  linear  ship  motion  theory  (along  with  strip  theory)  leads  to  the 
computation  of  the  response  amplitude  operator  as  a  function  of  relative 
wave  incident  angle,  frequency,  and  ship  speed.  Then,  under  the  assumption 
that  the  theory  of  linear  superposition  over  the  frequency  domain  holds  true, 
this  function  can  then  be  used  as  the  "black  box"  by  which  the  output 
spectrum  is  obtained  from  the  input  (wave  energy)  spectrum. 

(cD, )  =  |£/\ (C0„  Po ,  )r  •  (to. .  M 0 .  ) ; 

where:  p  =  the  relative  angle  between  the  ship's  forward  motion  and  the 
dominant  icident  wave  direction; 

Uq  =  the  vessel's  forward  speed; 


cOg  =  the  wave  encounter  frequency  = 


CO- 

0) 

•cospo 

8  ) 

(to^jPo,!/^) 


S,(co) 

l-(2o)(/o/g)-cospo| 


In  this  study,  information  on  relative  heading  anlgle  is  not  available; 
therefore,  it  will  be  assiuned  that  the  "worst  case"  relative  heading  in  view  of 
vertical  wave  bending  moment  corresponds  to  either  direct  head  or  following 
seas  (po=  0“,  180°). 


The  calculation  of  extreme  values  of  VWBM  which  will  be  discussed  in 
section  3.5.5  requires  that  the  value  of  the  average  (or  expected)  vertical 
wave  bending  moment  for  a  specific  sea  condition  be  known.  From  spectral 
analysis,  the  area  under  Svwbm(“)  or  the  zero^**  moment  of  Svwbm(o>)»  nk)  is 
equal  to  the  mean  square  value  of  the  response  (Evwbm)* 
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And  the  average  value,  VWBM ,  is  related  to  the  mean  square  value  by  the 
expression; 

VH™(//„r,)  =  0.866V£;:j^ 

Thus,  for  a  particular  vessel,  a  table  similar  table  can  be  developed, 
detailing  the  average  vertical  wave  bending  moment  response  to  a  given 
combination  of  Hg  and  T^.  This  table  would  represent,  for  a  specific  vessel,  a 
complete  set  of  input  data  for  the  purpose  of  calculating  the  extreme  values 
of  vertical  wave  bending  moments  that  the  vessel  might  encoxmter  during  its 
lifetime. 

3.4.4.  Stillwater  Bending  Moment 

It  should  be  remembered  that  a  tanker  typically  divides  a  significant  amount 
of  its  "at  sea"  time  between  at  least  two  different  loading  conditions.  In  this 
study,  both  full  load  and  ballast  conditions  will  be  treated  for  each  vessel. 
This  distinctions  affects  not  only  the  Stillwater  moment,  but  also  the  response 
of  the  vessel  to  wave  action  due  to  perhaps  a  different  draft  line  or  more 
significantly,  different  inertia  effects  resulting  fi*om  a  redistribution  of 
weight  from  one  loading  condition  to  the  other. 

The  calculation  of  SWBM  is  a  simple  matter  of  hydrostatics  and  involves  the 
difference  between  the  Weight  and  Buoyancy  distributions  along  the  length 
of  the  vessel.  Although  the  Stillwater  bending  moment  can  be  controlled  to  a 
certain  extent  and  calculated  fairly  accurately,  there  still  remains  a 
significant  element  of  imcertainty  in  its  representation  for  analytical 
purposes.  Nikolaidis  and  Kaplan  [5]  analyzed  data  presented  by  Guedes 
Soares  and  Moan  (1988)  and  predicted  that  the  standard  deviation  of 
Stillwater  bending  moments  for  a  particiilar  tanker  is  about  0.21  multiplied 
by  the  rule  based  value.  Treating  the  ballast  and  full  load  conditions 
separately  would  do  much  to  reduce  this  estimation  of  uncertainty. 
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Due  to  the  fact  that  the  tanker  operators  have  a  fair  amount  of  control  and 
information  about  the  Stillwater  bending  moment,  the  description  for  SWBM 
for  both  loading  conditions  will  be  left  up  to  the  user  to  supply  at  the 
beginning  of  the  life  assessment  rather  than  derived  from  fleetwide  data  as 
done  by  Nikolaidis  and  Kaplan.  It  will  be  assumed  that  SWBM  follows  a 
normal  (Gaussian)  probability  law  given  by: 


where  S  is  a  random  variable  representing  SWBM,  m  is  equal  to  the  mean 
value,  SWBM ,  and  Cg  is  the  standard  deviation.  Thus,  the  two  values,  m  and 
Gg,  need  only  be  supplied  in  order  to  describe  SWBM  for  a  particular  vessel. 

3.4.5.  Extreme  Total  Vertical  Bending  Moment  Distribution 

With  a  description  of  a  vessel's  environment,  response  to  the  environment  in 
terms  of  vertical  moment,  and  Stillwater  bending  moment  characteristics 
established  as  in  the  preceding,  the  extreme  value  distribution  of  the  Total 
Vertical  Bending  Moment  can  be  developed  for  both  full  load  and  ballast 
loading  conditions. 

The  basic  time  increment  involved  in  this  study  is  a  one-year  period.  This 
constitutes  a  "long-term"  situation  in  view  of  ocean  statistics.  While  this  fact 
does  not  affect  the  Stillwater  component  of  the  total  vertical  bending  moment, 
it  carries  strong  implications  for  the  interpretation  of  the  wave  statics  and 
vessel  response.  "Long-term"  implies  that  the  vertical  wave  bending 
response  of  a  vessel  during  this  time  period  can  not  be  described  by  a 
stationary  statistical  model.  However,  empirical  studies  have  shown  that  the 
amplitude  of  the  vertical  wave  bending  response  over  the  long  term  follows 
(approximately)  an  exponential  probability  law  with  the  average  (expected) 
value  of  the  wave  bending  moment  as  a  parameter: 
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x^O 

where  X  is  a  random  variable  representing  VWBM,  and  X,  =  VWBM  = 
average  value  of  VWBM. 

Defining  the  random  variable,  as  the  extreme  value  of  wave  bending 
moment,  X,  in  n  records  of  X,  the  use  of  order  statistics  permits  the 
probability  law  which  governs  to  be  expressed  as  follows  [Mansour,  JSR 
72]: 


As  stated  before,  Hogben  and  Lumb's  data  were  collected  over  a  period  of 
about  seven  years  and  therefore  constitutes  a  seven  year  record.  The 
parameter  n  in  the  above  equation  can  be  estimated  for  a  particular  vessel  as 
the  nearest  integer  to  the  value  of  the  vessel’s  design  life  (in  years)  divided 
by  seven.  For  example,  most  vessels  have  a  design  life  of  approximately 
twenty  years  and  consequently,  they  span  roughly  three  record  periods  of 
Hogben  and  Lumb's  sea  data;  i.e.  n  «  3. 

Therefore,  in  order  to  completely  know  the  distribution  of  the  extreme  value 
of  vertical  wave  bending  moment,  the  value  of  the  average  wave  bending 
moment  over  the  seven  year  record  period  is  the  only  remaining  item  to  be 
calculated.  Given  that  the  average  response  to  each  sea  state  has  been 
calculated  along  with  the  probability  that  the  vessel  will  experience  that  sea 
state,  the  total  average  wave  bending  moment  is  then  simply: 

H.  T, 
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where  y  corresponds  to  either  Ballast  (1)  or  Full  Load  (2)  conditions 

Letting  T  be  a  random  variable  representing  the  total  vertical  bending 
moment,  the  equation  for  TVBM  can  now  be  expressed  as: 

T„  =  S±Yn 

Combining  the  two  probability  laws  governing  S  and  Yn,  leads  to  the 
following  expression  for  T  (see  Mansour,  1972  for  derivation): 


dt 


— i— 


1-e 


-in-l 


Theoretically,  this  process  would  be  repeated  for  each  section  of  the  vessel  in 
both  loading  conditions  which,  for  the  case  of  a  vessel  with  twenty  designated 
stations,  would  result  in  forty  repetions  of  a  process  that  is  already 
computationally  demanding.  In  order  to  reduce  this  demand,  it  is  possible  to 
develop  the  above  expression  for  just  the  midship  section  in  each  of  the 
loading  conditions  and  then  make  assumptions  as  to  how  TVBM  varies  along 
the  length  of  the  vessel.  For  example,  both  the  mean  of  TVBM  and  the 
variance  could  be  assumed  to  vary  along  a  vessels  length  according  to  a 
distribution  factor  illustrated  in  Fig.  3.12: 


Disuibution 

Factor 


Figure  3.12:  Distribution  factor  of  TVBM  along  the  vessel  length 
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While  this  may  be  a  bold  assumption,  it  reduces  the  necessary  repetitions 
from  forty  to  two  (one  for  each  load  case). 

3.4.6.  Local  Loads 

The  calculation  of  local  loads  (i.e.  axial  stresses  on  each  element  as  defined 
for  a  section’s  structure)  can  be  acheived  by  employing  beam  theory.  The 
axial  stress  on  an  element  at  station  x,  and  at  a  distance  fi*om  the 
instantaneous  neutral  axis  is  given  by: 


TVBM{xyy„, 

I{x) 


where  I(x)  is  equal  to  the  area  moment  of  inertia  of  the  section  about  the 
neutral  axis. 
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3.5.  Capacity  Module 


Purpose:  To  produce,  for  a  specified  vessel,  probability  distribution  functions  of 

the  capacity  of  the  vessel’s  structure  to  resist  the  failure  inodes  defined  in 
the  Failure  Definition  Module. 

3.5.1.  Overview 

As  the  second  aspect  to  the  reliabiUty  problem,  the  Capacity  module 
generates  a  probabilistic  description  of  a  vessel’s  structure  to  resist  the 
seaway  loading  in  both  the  hogging  and  sagging  modes.  The  capacity  of  the 
structures  defined  in  the  Vessel  Description  module  can  be  generally 
described  m  terms  of  their  load/displacement  curves.  This  applies  to  both 
local  and  ultimate  failure  modes. 
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Material 

Stress-Strain 


2.1.6.1 


The  event  of  an  element  surpassing  its  elastic  limit  due  to  extreme  loading 
does  not  necessarily  constitute  failure  as  there  is  generally  some  residual 
strength  left  in  the  plastic  regime.  In  order  to  maintain  generality  in  this 
study,  structural  capacity  will  be  defined  as  that  level  of  load  at  which  the 
slope  of  the  load/displacement  curve  reaches  zero.  This,  in  effect,  defines  the 
ultimate  limit  state  for  each  element  and  group  of  elements.  While  some 
elements  may  need  to  be  replaced  or  repaired  at  lower  limit  states  such  as 
the  elastic  limit  or  some  limiting  value  on  displacement,  this  information  is 
too  detailed  for  the  general  treatment  presented  here. 

What  follows  is  the  development  of  these  load/displacement  curves  for  the 
structures  defined  in  the  Vessel  Description  Module. 

3.5.2.  Element  Load/End-Shortening  Curves 

Specifically,  the  load/displacement  curve  for  a  particular  element  is  cast  in 
terms  of  axial  load  vs.  the  shortening  at  the  ends  of  the  element.  In  view  of 
buckling,  there  is  a  high  degree  of  geometric  non-linearity  involved  in  the 
computation  of  this  relationship.  While  there  are  many  design  equations  and 
theories  available  to  predict  these  curves  and  the  buckling  capacity  of 
stiffened  panels,  these  generally  have  as  their  basis  a  linear  formulation  with 
some  correction  factor  to  account  for  non-linearities. 

The  most  rational  approach  is  to  deal  with  the  non-linearities  directly  in  a 
non-linear  finite  element  formulation.  With  the  proper  load  and  boundary 
information  supplied  by  the  Vessel  Description  module,  the  structural 
response  of  each  element  can  be  determined. 

Shown  below  are  some  possible  finite  element  models  that  can  be  used  to 
describe  the  types  of  structural  arrangements  and  response  behavior. 
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plate  elements 


plate  elements 

beam  elements 


orthotropic  plate  elements 

u'  o^T^  rir" 


Figure  3.15:  Possible  Finite  Element  Models 
There  are  four  types  of  general  responses  that  need  to  be  modeled.  They  are: 

•  buckling  of  plating  between  stiffeners 

•  column  buckling  (stiffener  or  plate  induced) 

•  stiffener  tripping  (or  torsional  buckling) 

•  overall  panel  buckling 


For  a  given  element  geometry,  each  response  mode  might  require  a  distinct 
finite  element  model  in  order  to  accurately  reproduce  the  intended  structural 
behavior.  This  could  potentially  result  in  a  large  computational  effort, 
especially  if  the  number  of  different  types  of  elements  for  the  section  under 
consideration  is  large. 
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As  an  alternative,  the  design  equations  mentioned  previously  can  be  used, 
combining  both  analytical  theory  and  empirical  data  to  predict  the  critical 
stresses  at  which  a  stiflfened  panel  might  buckle  and,  further,  to  develop 
approximate  load/end-shortening  curves  for  a  given  element  cross  section. 
What  follows  is  an  illustrative  procedure  for  determining  these  curves  based 
on  design  equations. 

It  can  be  assumed  that  the  load/end-shortening  curve  for  a  particular 
element  can  be  based  on  the  stress-strain  curve  of  the  material  (steel)  of 
which  it  is  made.  The  element  stress-strain  curve  will  follow  the  material 
stress-strain  curve  up  until  a  critical  point  at  which  the  element  becomes 
unstable  in  compression  and  buckles.  It  is  necessary,  therefore  to  first 
develop  material  stress-strain  curves. 

Generalized  material  stress-strain  curves  can  be  developed  based  on  a 
relatively  small  number  of  parameters.  More  specifically,  given  (for  a 
particular  steel)  the  elastic  section  modulus,  E,  the  yield  stress,  Oyp,  the 
proportional  limit  stress  Op,  and  Poisson's  ratio  v,  it  is  possible  to  estimate 
the  stress  vs.  strain  characteristics  of  that  steel  in  a  complete  yet 
approximate  sense. 

This  is  achieved  by  dividing  the  stress-strain  curves  into  three  regions 
signifying:  (I)  the  linearly  elastic  range,  (II)  the  nonlinear  elastic  range,  and 
(III)  the  perfectly  plastic  (yield)  range  as  shown  in  the  figure  below. 


Steel  Type 

E 

Elastic  Modulus 

Yield  stress 

% 

Proportional  limit 

V 

Poisson's  ratio 
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Figure  3.16:  Material  Properties 


Region  I  e<E^ 

The  stress-strain  relationship  in  this  region  is  simply  linear  with  the  Elastic 
modulus  as  the  constant  of  proportionality,  i.e. 

o  =  £e 

Region  II 

This  region  generally  signifies  a  nonlinear  "softening"  of  the  material  which 
is  represented  by  a  gradual  change  in  the  slope  of  the  stress-strain  curve. 
An  expression  for  the  curve  in  this  region  is  as  follows: 


/■‘Kl 

G"  -O  O 

_  ^yp  ^yp^p 

‘  £c 

yp 

K'f  — 

Region  III 

Perfect  plasticity  is  assumed  in  this  region  which  is  to  say  that  the  material 
can  no  longer  continue  to  support  increasing  load  and  deforms  plastically 

With  the  material  behavior  defined  as  above,  it  remains  to  determine  the 
critical  stress  at  which  an  element  becomes  imstable  and  buckles.  This 
critical  stress  can  be  assumed  to  divide  the  load  end  shortening  curve  into 
two  regions,  one  region  in  which  the  element  behaves  according  to  the  stable 
material  behavior,  and  the  other  in  which  the  element  rapidly  "sheds"  its 
load  indicated  by  a  negative  slope  in  the  load/end-shortening  curve.  In 
general,  different  elements  will  exhibit  load  shedding  characteristics  to 
varying  degrees,  but  this  phenomenon  is  very  difficult  to  formulate 
mathematically. 

Buckling  of  the  plate  between  stiffeners  does  not  necessarily  result  in  the 
failure  of  the  stiffened  panel.  However,  the  buckling  strength  of  the 
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stiffener/attached  plating  combination  is  strongly  affected  by  the  stiffeness  of 
the  plate  between  stiffeners,  and  buckling  in  this  region  can  lead  to  a 
sigmhcant  reduction  in  the  stability  of  the  column  type  stiffener/plating 

combination.  This  effect  can  be  modeled  in  terms  of  the  "effective  width" 
concept. 

Shown  in  Fig.  3.17  is  a  flow  chart  describing  the  calculation  of  the  critical 
buckling  stress  of  an  element.  This  model  includes  the  effect  of  buckling  of 
the  plate  between  stiffeners  by  considering  that  if  the  critical  stress  of  the 
plate  between  stiffeners  is  less  than  the  that  of  the  stiffener  and  attached 
plate  as  a  column,  then  only  an  "effective"  width  of  the  plate  should  be 
considered  in  the  computation  of  the  colunm  strength.  The  effective  width  is 
calculated  using  the  computed  critical  column  stress.  The  column  stress  is 
then  recalculated  using  the  new  width  of  the  attached  plating  and  an 
iterative  process  is  begun. 
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Strain 


Figure  3.17;  Flow  Chart  for  Element  Capacity  CalruInHnn 

Another  case  might  arise  in  which  the  critical  stress  of  the  element 
considered  as  a  column  is  greater  than  the  ultimate  stress  of  the  material  of 
which  it  is  composed.  The  element  may  then  be  considered  as  a  "hard  spot" 
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meaning  that  the  behavior  of  the  element  in  compression  follows  that  of  the 
material. 


3.5.3.  Ultimate  Capacity 

While  the  finite  element  treatment  of  individual  elements  and  specific 
element  groups  provide  accurate  predictions  of  their  load/end-shortening 
curves,  it  is  not  feasible  to  apply  these  same  methods  to  the  entire  section's 
structural  system.  The  computational  complexity  of  such  a  problem 
precludes  the  development  of  an  interactive  PC-based  computer  application. 

Instead,  the  ultimate  capacity  of  a  given  section  can  be  determined  by  a 
method  that  was  outlined  by  Smith  [6].  This  general  procedure,  outlined 
below,  determines  a  Resisting  Moment  vs.  Curvature  relation  for  a  prismatic 
box  girder  section  based  on  the  individual  load/end-shortening  curves  of  its 
constitutive  structural  elements. 

Step  1;  determine  the  properties  of  the  section  (as  built  or  corroded).  This 
includes  calculating  the  second  moment  of  area,  I,  and  the  position  of  the 
neutral  axis. 

Step  2:  determine  the  elastic  limit  of  the  section,  i.e.  the  moment  and 
corresponding  curvature  at  which  the  first  element  in  the  cross  section 
reaches  its  elastic  limit  defined  by  its  individual  load/end-shortening  curve. 
Theoretically,  the  moment  curvature  relation  of  the  entire  section  is  linear  to 
this  point. 

Step  3:  firom  the  elastic  limit  point,  apply  an  incremental  curvature  to  the 
section  the  magnitude  of  which  can  be  arbitrarily  defined  as  a  small 
percentage  of  the  elastic  limit  curvature. 

Step  4;  with  the  assumtion  that  plane  sections  remain  plain,  the  strain  on 
each  element  can  be  calculated  as: 

•  C  ;  =  distance  of  element  above  neutral  axis 

C  =  section  curvature 
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Step  5:  from  each  element's  load/end-shortening  curve,  the  element  forces  can 
be  calculated  corresponding  to  the  strain  calculated  in  step  4. 

Step  6:  Since,  at  this  point,  at  least  one  element  has  passed  its  linear  elastic 
point,  a  softer  more  flexible  local  structure  will  result  in  an  imbalance  in 
horizontal  forces  and  a  shift  in  the  neutral  axis  is  required  to  ensure  that 
only  a  pure  bending  moment  is  acting  on  the  section.  This  generally  would 
require  an  iterative  procedure  where  the  neutral  axis  is  shifted  away  from 
the  plastic  region,  the  strains  and  forces  are  recalculated,  and  the  process  is 
repeated  until  there  is  a'  zero  net  horizontal  force.  However,  if  the  section 
curvature  increment  is  small  enough,  one  incremental  shift  of  the  neutral 
axis  can  be  assumed  to  be  accurate  enough,  and  is  given  by: 

SHIFT  = 

Cl(£,A,) 

Step  7.  determine  the  bending  moment  that  corresponds  to  the  current  state 
of  curvature.  This  moment  is  simply  calculated  as: 

A/  =  y  a  -4  -v 

c  €  e 

Step  8:  apply  the  next  increment  in  section  curvature. 

From  this  point  steps  4  through  8  are  repeated  xmtil  the  complete 
moment/curvature  relation  is  obtained.  The  flow  chart  in  Fig.  3.18 
graphically  illustrates  the  procedure. 
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Curvature 


3.6.  Corrosion  Module  (Time  Variability  of  Capacity) 


Purpose:  To  provide  a  means  by  which  to  calculate  corrosion  rates  from  general 
wastage  data  drawn  from  the  entire  fleet,  to  assign  these  corrosion  rates 
to  the  elements  and  sub>elements  of  the  specific  vessel  in  question,  and 
then  to  control  the  time  step  procedure  in  the  life  assessment. 

Although  the  inclusion  of  corrosion  rates  and  the  effect  of  corrosion  in  the  life 
assessment  procedure  is  a  fairly  straightforward  matter,  it  is  of  extreme 
importance  in  that  it  constitutes  the  time  variability  component  without 
which  there  would  be  no  life  assessment;  The  corrosion  module  consists  of 
three  parts: 


•  The  collection  of  corrosion  data, 

•  The  statistical  analysis  of  corrosion  data,  and 

•  The  integration  of  the  results  into  the  life  assessment  procedure. 

3.6.1.  Corrosion  Data  Collection  and  Modelling 

There  is  an  abundance  of  gauging  reports  from  which  data  can  be  drawn, 
collected  during  regular  inspections  of  the  entire  tanker  fleet  over  many 
years.  The  challenge  involved  in  this  part  of  the  Corrosion  Module  is  how  to 
model  the  data  in  such  a  way  that  trends  can  be  identified  that  will  be  useful 
to  the  analysis.  It  is  not  sufficient  or  rational  to  provide  just  one  number  as  a 
representation  of  the  corrosion  rate  situation  for  an  entire  vessel.  There  are 
many  factors  that  influence  the  wastage  of  tanker  structures  and  the  values 
for  mean  rate  can  vary  substantially  throughout  the  body  of  a  vessel.  Pollard 
[8]  compiled  the  following  list  of  important  factors  effecting  corrosion  rates: 


Ship  size 
Delivery  date 
Cargo  type 
Double  bottom 
Double  side 
Class  society 
Trade  route 
Tank  location 


Tank  type 
Time  in  cargo 
Time  in  ballast 
Corrosion  protection 
system 
Ballast  type 
Tank  temperature 
Tank  humidity 
Inert  gas 


Cargo  sulphur  content 
Cargo  water 
Wax  in  cargo 
Heated  cargo 
Tank  washing 
Corrosion  type 
Corroded  detail 
Location 
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The  corrosion  rates  to  be  used  in  the  life  assessment  procedure  are 
determined  through  a  statistical  analysis  of  the  corrosion  gauging  data 
stored  in  the  referenda!  database.  During  the  first  year  of  the  Structural 
Maintenance  for  New  and  Existing  Ships  Project  conducted  at  the  University 
of  California  at  Berkeley,  these  data  were  collected  and  anadyzed  [  ]. 
Corrosion  rates  were  categorized  by  a  combination  of  tank  type  and  detail 
type  and  also  by  a  combination  of  tank  type  and  general  location  within  the 
tank.  The  tank  types  that  are  considered  in  the  study  fall  into  the  following 
four  descriptions: 

1)  Cargo  only 

2)  Ballast  only 

3)  Cargo/clean  ballast 

4)  Cargo/dirty  ballast 

The  second  category  which  involves  trends  in  corrosion  wastage  as  they  are 
affected  by  general  location  within  the  tank  (upper  third,  middle  third,  lower 
third,  etc.)  provides  qualitative  information  only,  and  therefore  can  not  easily 
be  used  as  input  in  the  analysis.  While  location  within  the  tank  has  a 
significant  influence  on  the  corrosion  rate  of  the  structural  components,  the 
data  is  not  detailed  enough  to  provide  a  quantification  of  these  trends. 

A  further  deficiency  in  this  data  model  arises  when  one  considers  that  the 
tank-type/detail-type  category  only  gives  information  regarding  the  tank-type 
on  one  side  of  the  plating,  generally  the  side  on  which  the  longitudinal 
stiffeners  are  located.  An  area  of  longitudinal  bulkhead  plating,  for  example 
can  have  heated  cargo  on  one  side  and  cold  water  ballast  in  the  adjacent 
wing  tank,  while  another  area  of  longitudinal  bulkhead  plating  of  the  same 
cargo  tank  can  have  more  heated  cargo  in  the  wing  tank  on  the  other  side. 
This  situation  could  result  in  a  significantly  different  corrosion  rate  for  what 
would  be  considered  an  identical  tank-type/detail -type  combination  by  the 
database. 

In  developing  this  module,  efficient  use  can  be  made  of  the  way  in  which  the 
Vessel  Definition  Module  handled  the  input  of  each  structural  element.  A 
"key"  identifier  can  be  assigned  to  a  particular  sub-element  plating  at  the 
time  that  gauging  data  (thicknesses)  are  entered  into  the  database.  For 
example,  when  a  vessel  xmdegoes  inspection  and  a  measurement  is  taken  of 
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the  thickness  of  a  particular  sub-element,  knowledge  of  the  element  to  which 
the  plating  belongs  in  conjunction  with  the  section  number  corresponding  to 
that  element  will  give  access  to  the  following  information  from  the  element 
database; 


•  type  of  element  (side  shell,  longitudinal  bulkhead,  etc.) 

•  the  contents  of  the  tanks  to  either  side  of  the  plate  sub-element  (cargo  only, 
ballast  only,  etc.) 

•  the  region  within  the  tank  where  the  element  is  located  (ullage,  middle, 
lower,etc) 

•  the  original  thicknesses  of  the  sub-elements. 

This  information,  plus  the  identification  of  the  sub-element  type  (plate,  web, 
flange)  at  the  time  the  meastirement  is  input,  can  lead  to  a  fairly 
comprehensive  description  of  the  major  factors  that  are  involved  in  the 
identification  of  corrosion  rate  trends. 

3.7 .  Reliability  Module  (Probability  of  Failure) 

Purpose:  To  calculate  the  probability  of  failure  in  each  failure  mode  defined  in  the 
Failure  Definition  Module  based  on  the  demand  and  capacity  determined 
for  the  vessel  in  question. 

Three  major  sources  of  rmcertainty  in  the  failure  probability  calculations 
come  from  the  Capacity,  Demand,  and  Corrosion  modules.  In  this  study 
there  are  two  levels  of  structural  failure  to  be  examined;  the  element  (local) 
failure  level  and  the  ultimate  (global)  failure  level.  With  the  probability 
density  function  for  load  approximated  by  the  normal  distribution,  and  with 
the  corrosion  and  capacity  information  similarly  described,  the  entire 
reliability  problem  reduces  to  the  fundamental  level. 

For  a  particular  failiu-e  mode  i,  a  "safety  margin",  M  can  be  defined  as: 

The  probability,  then,  that  M<0  is  equal  to  the  probability  that  the  capacity 
of  the  structure  to  resist  failure  mode  i  is  less  than  the  demand  that  is  placed 
on  the  structure,  which  in  turn  is  simply  the  probability  of  failure  in  mode  i. 
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M  constitutes  a  random  variable,  also  normally  distributed,  whose  mean 
value,  and  standard  deviation,  o,  can  be  easily  calculated  for  any  time 
instant  and  any  failure  mode.  Assuming  independence  between  the  capacity 
and  the  demand: 


The  probability  of  failure  is  therefore: 

,  =  p[M  <  0]  =  (bf  = <b(-p) 

\  ) 

where  p  =  is  defined  as  the  "safety  index"  and  can  be  thought  of  as  the 
number  of  standard  deviations  by  which  exceeds  zero. 

3.8.  Availabiuty 

With  the  various  global  variables  defined  as  they  have  been  in  the  preceding, 
there  are  a  number  of  ways  in  which  the  reliability  of  a  particular  vessel  can 
be  formulated  within  the  context  of  its  life  span.  For  example,  the  reliability 
can  be  cast  in  terms  of  the  mean  and  standard  deviation  of  the  time  that  it 
takes  for  the  calculated  availability  to  drop  below  the  design  value,  or, 
alternatively,  the  mean  and  standard  deviation  of  the  availability  can  be 
presented  at  the  end  of  tiie  vessel's  design  life.  These  two  alternative 
formulations  are  presented  in  figure  3.20.  The  subtle  difference  between  the 
two  is  that  in  the  first  formulation,  the  uncertainty  in  the  time  dimension  is 
treated  while  the  limiting  availability  is  taken  as  deterministicCdesign  Av), 
whereas  in  the  second  formulation,  the  weight  of  uncertainty  rests  on  the 
availability  dimension  while  the  design  life  determines  the  limiting  time. 
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Av  iyn/yr) 


Av  iyn/yr) 


For  the  purpose  of  defining  corrosion  limits,  it  is  more  important  that  the 
uncertainty  in  time  is  treated  since  the  time  dimension  is  directly  involved  in 
the  determination  of  corrosion  rates. 
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4.  Synthesis  of  the  Model  Application 


4.1.  Overview 

Presented  in  this  chapter  is  the  initial  development  of  a  computer  application 
which  is  modelled  after  the  approach  detailed  in  the  preceding  chapters.  The 
following  simplifications  were  made  to  facilitate  the  development  of  the 
model: 


1)  The  load  definition  module  was  eliminated,  using  in  its  place  a  rule 
based  definition  of  the  extreme  bending  moment  amidship. 

2)  The  'time  until  failure’  was  based  solely  on  the  ultimate  capacity  of 
a  section  to  withstand  the  imposed  bending  moment.  While  this  is 
only  a  component  of  a  component  of  the  overall  availability  of  a  vessel, 
it  clearly  and  adequately  represents  a  limiting  condition,  and  further, 
the  additional  data  and  routines  needed  to  complete  the  entire 
availability  calculation  was  beyond  the  scope  of  this  project. 

3)  The  database  files  are  accessed  directly  through  FOXPRO, 
eliminating  the  need  for  the  development  of  input  screens. 

4)  Corrosion  rates  were  'hard-wired'  or  manually  input  to  the  system 
due  to  the  fact  that  the  data  did  not  exist  in  the  designed  format. 
Corrosion  rates  were  based  on  Pollard's  findings. 

Shown  in  Figure  4.1  is  the  context  layer  diagram  for  the  model  SMIS 
application.  When  compared  with  Figure  1.1  this  diagram  expresses  the 
above  simplifications  in  graphic  form.  By  employing  a  rule  based  definition 
of  the  loads  and  by  eliminating  the  need  for  all  Uiat  is  required  for  the 
calculation  of  the  three  components  of  Unavailability,  the  external  inputs 
become  simply  the  midship  section  idealization  and  the  general  pareimeters 
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used  in  the  calculation  of  the  extreme  midship  bending  moment.  The  system 
then  bases  its  definition  of  wastage  limits  on  the  prediction  of  the  time  that  it 
takes  for  the  wastage  of  the  scantlings  to  decrease  the  ultimate  capacity  of 
the  section  to  the  point  at  which  it  no  longer  can  satisfactorily  withstand  the 
rule  based  load. 


Figure  4.1:  Context  layer  diagram  for  the  model  SMIS  application 


4.2.  Database  Management 

Without  the  need  to  support  the  entire  unavailability  calculation,  the 
structuring  and  management  of  the  database  becomes  considerably  simpler. 
Direct  input  of  data  to  the  database  files  further  reduces  the  complexity  of 
the  database  management  issue  by  eliminating  the  need  for  input  screens. 
The  principal  components  of  the  database  are  simply  the  CLASS  and 
SECTION  idealization  databases.  Shown  in  Figure  4.2  is  the  next  layer  DFD 
followed  by  the  structure  of  the  developed  database  (Fig.  4.3). 
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Fig  4.2:  Second  Laver  DFD 
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application 


The  definition  of  a  vessel  is  simply  a  matter  of  filling  the  data  structure 
defined  above. 


4.3.  The  Analytical  Session 

4.3.1.  Setting  up  an  Analytical  (Life  Assessment)  Session 

Setting  up  an  analytical  session  involves  little  more  than  identifying  which 
section  of  which  vessel  is  to  be  considered  and  any  biases  that  are  applied  to 
customize  the  loading  condition.  In  addition,  the  time  increment  in  years 
needs  to  be  designated  along  with  a  minimum  value  for  the  safety  index,  p 
upon  \\  hich  the  life  assessment  is  based.  A  single  main  screen  was 
developed  to  accept  the  input  of  these  session  parameters,  do  some 
elementary  calculations  and  prompt  the  user  to  begin  the  analysis  once  the 
session  parameters  have  been  defined.  Shown  in  Fig.  4.4  is  this  Main  Screen 
as  it  appears  on  the  monitor. 
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Values  for  the  Stillwater  and  vertical  wave  bending  moments  in  both  the 
hogging  and  sagging  condition  are  automatically  computed  once  a  particular 
vessel  class  is  chosen  from  the  popup.  At  the  same  time,  the  "sections"  popup 
is  filled  with  section  names  for  the  chosen  vessel  and  the  input  fields  for 
biases  and  coefficients  of  variation  are  enabled  and  defa\ilt  values  are 
displayed  (1.0  for  biases,  and  0.0  for  COV’s).  A  field  also  exists  for  a  user 
input  "session  id  code"  which  is  stored  in  the  support  database  and  can  be 
used  to  distinguish  the  results  of  a  particular  anal3rtical  session  for  later 
study. 

Once  all  of  this  information  has  been  input,  including  the  time  step  (typical 
values  should  be  aroxmd  five  to  seven  years  for  this  increment)  and  minimum 
p  level,  a  button  labeled  "Begin’  starts  the  analysis. 

4.3.2.  Analysis 

The  main  analysis  routines  lie  beneath  the  setup  screen  and  automatically 
control  the  flow  of  the  analysis,  the  links  to  the  support  database,  and  the 
generation  of  results.  The  program  calculates  the  safety  index  at  a  particular 
time  step  based  on  the  computation  of  the  designated  section’s  ultimate 
capacity  in  both  hogging  and  sagging  modes  and  the  combination  of  the 
defined  Stillwater  and  vertical  wave  bending  moments.  If  the  calcvdated 
safety  index  is  greater  than  the  defined  minimum,  then  the  "age"  of  the 
vessel  is  increased  by  one  time  increment,  the  section’s  scantlings  are 
reduced  according  to  the  appropriate  corrosion  rates,  and  the  process  is 
repeated  until  the  safety  index  drops  below  the  set  minimum.  In  addition  to 
calculating  the  safety  index  at  each  time  step,  the  program  builds  moment 
curvature  diagrams  according  to  the  procedure  outlined  in  chapter  three. 

The  following  is  a  list  of  the  principle  modules  that  comprise  the  analysis 
routine  and  a  brief  description  of  each. 

Module: _ MAIN 

Purpose:  This  is  the  main  module  that  performs  the  remaining 

preliminary  computations  regarding  loads,  coordinates  the 
subroutines,  and  generally  controls  the  flow. 

Input:  Session  control  parameters  from  main  screen 
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Output:  Global  information  regarding  the  section  at  each  time  step 

which  is  stored  in  the  LIFE.DBF  database  (i.e.  safety  index,  hog 
capacity,  sag  capacity,  initial  neutral  axis,  etc) 

Called  by:  Main  Screen 


Calls: 


CORRODER,  CAPACITY 


Purpose: 


Input: 


Output: 


This  module  develops  the  important  "elms"  array  which  contains 
the  element  specific  information  such  as  dimensions  and 
material  properties  and  includes  the  correct  plate  thicknesses 
according  to  current  age  of  teh  vessel  and  the  appropriate 
corrosion  rates. 

"Temp"  array  which  is  downloaded  from  the  ELEMENTS.DBF 
database  and  contains  the  element-by  element  description  of  the 
chosen  section. 

Elms  array  which  is  similar  to  the  Temp  array  but  has 
updated  the  element  dimensions  to  reflect  wastage  of  an  aged 
vessel. 


Called  by:  MAIN 


Calls: 


none 


Purpose:  This  module  follows  the  procedure  for  calculating  ultimate 

capacity  of  a  section  based  on  the  load/end-shortening  curves  of 
its  constituant  elements.  Program  is  excecuted  for  a  particular 
time  step. 

Input:  "Elms"  array  described  above 

Output:  caphog  and  capsag  (ultimate  capacity  in  hogging  and  sagging 

conditions)  both  of  which  are  returned  to  the  MAIN. 
Additionally,  the  moment  curvature  relationship  is  uploaded  to 
MOMCURV.DBF. 

Called  by:  MAIN 
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Calls; 

ELMPROPS,  STRSSTRN 

Module: _ 

ELMPROPS 

Purpose; 

This  module  calculates  three  element  properties  which  are 
added  to  the  elms  array. 

Input: 

The  information  contained  in  one  record  (row)  from  the  elms 
array 

Output: 

crjstrs  (critical  buckling  stress),  area  (cross  sectional  area),  and 
inertia  (moment  of  inertia  about  the  centroid)  particular  to  an 
element,  and  stored  in  columns  14,15,  and  16  of  the  elms  array 

Called  by: 

CAPACITY 

Calls: 

none 

Module: _ 

Purpose: 

This  module  returns  the  stress  corresponding  to  an  input  strain 
level  for  a  particular  element  based  on  its  material  properties 
and  critical  buckling  stress. 

Input: 

crjstrs,  elasmod,  ujstrs,  yjstrs,  poisson,  strain 

Output: 

stress 

Called  by: 

CAPACITY 

Shown  in  Fig.  4.5  is  a  schematic  view  of  the  modules  involved  in  the  initial 
application. 
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Figure  4.5:General  Flow  Chart  for  the  Analysis 

While  the  CAPACITY  module  is  fairly  robust  and  genuine  in  its  approach, 
the  ELMPROPS  module  which  calculates  the  capacity  of  an  element  is 
scarcely  more  than  symbolic  and  the  development  of  a  sophisticated 
load/endshortening  module  would  add  significantly  to  the  value  of  this 
program. 
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5.  Example  Application 


5.1 .  Example  Vessel 


The  vessel  that  was  chosen  for  the  example  application  is  a  216,000  DWT 
single  bottom  VLCC  named  the  Energy  Concentration.  In  July  of  1980,  the 
Energy  Concentration  suffered  a  "broken  back"  while  discharging  oil  at  the 
Mobil  Terminal  in  Rotterdam.  While  there  were  many  factors  that  lead  up  to 
the  ultimate  collapse,  the  fact  that  the  VLCC  was  ten  years  old  at  the  time 
suggests  that  wastage  of  the  structure,  particularly  the  bottom  plating  and 
longitudinals,  must  have  played  an  important  role.  In  addition  to  presenting 
an  interesting  corrosion  study,  this  event  was  extensively  studied  by 
Rutherford  and  Caldwell  [Ultimate  Longitudinal  Strength  of  Ships:  a  Case 
Study]  the  results  of  which  can  be  used  as  a  comparison. 

A  brief  description  of  the  physical  characteristics  of  the  Energy 
Concentration  is  given  in  the  following  tables  and  figures. 


General  Particulars 


L.O.A. 

L.B.P. 

Breadth  (mid) 
Depth  (mid) 
Gross  tonnage 
Deadweight 
Block  Coef. 


326.75  m 
313.0  m 
48.19  m 
25.2  m 
98,894  tons 
216,269  tons 
0.809 


The  overall  design  and  layout  of  the  Concentration  was  typical  of  VLCC's 
built  around  1970.  The  cargo  section  of  the  hull  was  divided  by  two 


70 


longitudinal  and  seven  transverse  bulkheads,  making  a  total  of  five  center 
tanks  and  twelve  wing  tanks  (Fig.  5.1). 


Figure  5.1;General  Arrangement  (Profile  and  Plan) 


The  catastrophic  collapse  of  the  Concentration  occured  around  frame  76 
where  the  Stillwater  bending  moment  was  at  its  maximum  value  of  roughly 
17,940  MNm.  Shown  in  Fig.  5.2  is  a  section  view  of  the  longitudinal 
structure  of  this  mid  body  portion  of  the  hull.  The  bottom,  side,  deck  and 
longitudinal  bulkhead  plating  are  reinforced  by  closely  spaced  longitudinal 
stiffeners.  The  logitudinal  structure  is  then  supported  by  transverse  web 
frames  spaced  5.1  meteres  apart. 
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For  this  study,  the  hull  was  idealized  as  a  combination  of  over  two  hundred 
plate  stiffener  element  combinations.  A  full  list  of  the  elem tents  that 
constitute  this  section  is  given  in  the  appendix  to  this  report.  In  the  table 
below  are  a  few  examples  of  the  input  required  to  define  an  element  of  the 
cross  section. 


Element  id: 

BOIS 

S17S 

D05S 

Configuration: 

H 

Plate 

plate  breadth  (mm) 

plate  thickness 

(mm) 

corrosion 

rate(mm/yr) 

1000 

25 

0.197 

925 

23.5 

0.051 

1000 

25.0 

0.11 

Web 

web  depth  (mm) 

797 

747 

web  thickness  (mm) 

15 

12.7 

32 

corrosion 

rate(mm/yr) 

0.063 

0.035 

Flange 

flange  width  (mm) 

200 

180 

- 

flange  thickness 

(mm) 

33 

25 

• 

corrosion 

rate(mm/yr) 

0.053 

0.050 

• 

Materia 

1 

Type 

HTS 

MS 

HTS 

elastic  mod 
(N/mm'^2) 

4233 

4233 

4233 

yield  stress 
(N/mm'^2) 

350 

315 

350 

ultimate 

strs(N/mm'^2) 

555 

525 

555 
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5.2.  Load  Condition 


The  Concentration  failed  while  in  port  and  in  the  hogging  condition.  The 
demand  placed  on  the  vessel  consisted  solely  of  a  hogging  still  water  bending 
moment.  In  order  to  simulate  this  demand  situation  for  the  SMIS 
calculation,  the  sagging  loads  can  be  left  as  they  were  calculated  since  these 
represent  non-extreme  loads,  and  are  unlikely  to  drive  the  overall  safety 
index.  The  hogging  vertical  wave  bending  moment  can  be  eliminated  by 
setting  its  bias  factor  equal  to  zero.  Finally,  the  extreme  Stillwater  hog 
moment  of  17,940  MNm  be  derived  from  the  "rule  based"  by  assigning  an 
appropriate  bias  factor.  The  Stillwater  hog  moment  calculated  by  the  SMIS 
was  5,851  MNm  implying  a  required  bias  factor  of  around  3.0. 
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6.  Results  and  Conclusions 


6.1 .  Results  of  the  Example  Appucation 

The  ultimate  capacity  of  the  Concentration  in  the  hogging  mode  is  shown  in 
Fig.  6.1.  The  capacity  was  calculated  for  each  of  six  time  steps  ranging  from 
the  zero  year  (as-built)  section  up  until  the  ten  year  mark  which  corresponds 
to  the  age  of  the  vessel  when  it  sailed  into  port  for  the  last  time.  The 
horizontal  line  in  each  graph  represents  the  extreme  Stillwater  load  applied 
to  the  reliability  calculation. 


Bending  Moment  MNm  Bending  Moment  MNm 


76 


847 


The  capacity  of  the  midship  section  structure  in  the  hogging  mode  was  found 
to  decline  almost  linearly  with  time  at  a  rate  of  roughly  1,500  MNm/year. 
This  trend  is  displayed  graphically  in  Fig.  6.2 


Figure  6.2;  The  Decline  in  Ultimate  Capacity  of  the  Midship  Section 


6.2.  Reliabiuty  of  Sample  Vessel  and  Corrosion  Limits 

Presented  below  is  the  decline  in  ’Tseta"  over  the  lifespan  of  the  Energy 
Concentration  as  calculated  by  in  the  example  application.  The  trend  that 
the  numbers  display  illustrates  what  would  be  expected.  The  graph 
represents  a  slightly  accelerated  (i.e.  non-linear)  decrease  in  the  safety  index. 
This  is  the  result  of  a  higher  percentage  of  elements  entering  the  non-linear 
regions  of  their  load/end-shortening  curves.  The  non-linearity  would  be  more 
pronounced  if  the  uncertainty  in  corrosion  rates  were  included  in  the  model 
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resulting  in  a  "spreading  out"  of  the  probability  density  function  for  the 
ultimate  capacity. 


Figure  5.4;  Decline  in  Reliability  over  Time 

6.3.  Corrosion  Limits 

While  there  where  many  factors  involved  in  the  failure  of  the  Energy 
Concentration,  clearly  the  wastage  of  the  internal  structure  was  one  of  them. 
The  wealth  of  information  surrounding  this  particular  event  provided  the 
insight  in  this  case  that  otherwise  would  have  to  come  from  a  very  thourough 
treatment  as  per  the  approach  outlined  in  this  study.  Only  a  comprehensive 
treament  of  all  the  major  factors  involved  will  lead  to  acciirate  predictions  of 
the  allowable  wastage  limits. 

In  the  mean  time,  with  the  benefit  of  hindsight  and  a  historically  based 
estimate  of  corrosion  rates,  wastage  limits  can  be  assigned  to  each  sub¬ 
element  of  the  failed  section.  Applying  a  ten  year  time-until-friilure,  and 
assuming  that  corrosion  rates  remain  constant  over  tlie  long  term,  tlie 
allowable  wastage  can  be  calculated.  The  results  for  the  bottom  shell 
element  presented  earlier  are  as  follows: 
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Subelement 

corrosion  rate 
(mm/yr) 

time-imtil  failure 
(yrs) 

allowable 
wastage  (mm) 

Plate 

0.197 

X 

10 

1.97 

Web 

0.063 

X 

10 

0.63 

Flange 

0.053 

X 

10 

0.53 

6.4.  Recommendations  for  Future  Work 

The  goal  of  this  project  was  to  develop  a  rational  approach  to  the  definition  of 
corrosion  limits  in  tankers.  The  basic  framework  for  an  analytical  tool  that 
can  be  used  to  solve  this  problem  has  been  laid  out  in  the  preceding  sections 
along  with  an  example  application  to  illustrate  the  procedure.  There 
remains  much  work  to  be  done  in  order  to  implement  the  ideas  presented  in 
this  study  and  create  an  application  that  can  be  used  by  the  industry. 

A  endeavor  such  as  this  one  actually  involves  two  disciplines:  Naval 
Architecture  and  Computer  Systems  Analysis.  In  order  to  develop  a  working 
application  based  on  this  study,  the  expertise  of  both  fields  are  required. 
Listed  below  are  reccomendations  for  future  work  in  both  areas: 

6.4.1.  Naval  Architecture  Topics 

•  Element  Behavior:  Among  the  many  uncertainties  involved  in  the 
modeling  of  this  problem,  the  prediction  of  the  individual  element 
behavior  plays  a  very  significant  role.  This  xmcertainty  alone,  if  not 
treated  properly,  could  potentially  invalidate  the  results  of  the  reliability 
calculation.  Topics  to  be  considered  are:  lateral  hydrostatic  pressure, 
initial  imperfections,  and  residual  stresses. 

.  RnnnHarv  Conditions  and  Interaction  Between  Elements:  The  effect  that 
adjacent  elements  have  on  each  other  plays  a  significant  role  in  the 
calculation  of  buckling  strength. 

•  T,nading  Module:  A  loading  module  based  on  ship  motion  theory  needs  to 
be  incorporated  into  the  application. 
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•  Uncertainty  and  Correlation:  A  comprehensive  treatment,  module  by 
module,  of  all  the  uncertainties  involved  is  crucial  to  the  success  and 
usefullness  of  the  system.  This  includes  modellini^  correlation  between 
failure  modes  and  correlation  between  repair  times. 

•  HQE:  Some  recent  work  being  done  in  the  area  of  Human  and 
Organizational  Errors  could  be  incorporated  into  the  approach. 

6.4.2.  Computer  Systems  Topics 

•  Data  Gathering:  In  order  to  support  all  of  the  aspects  of  the  life 
assessment  procedure,  the  database  must  have  sufficient  and  accurate 
data. 

•  Data  MQ.dgling:  This  task  involves  taking  the  raw  data  and  setting  it  in  a 
format  that  can  be  used  by  the  analytical  routines  while  at  the  same  time 
providing  for  the  ease  of  input. 

•  User  Interface:  A  consitent  user  interface  needs  to  be  designed  that  will 
allow  a  range  of  users  to  operate.  This  includes  providing  help  screens, 
menu  bars,  error  checking,  input  (data  entry)  screens,  and  output 
(reporting)  screens. 
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Appendix  A:  Source  Code  for  SMIS 


11/30/93 


ANALYSIS.SPR 


22:35:45 


Author's  Name 

Copyright  (c)  1993  Company  Name 

Address 

City,  Zip 

Description: 

This  program  was  automatically  generated  by  GENSCRN. 


#REGION  0 

REGIONAL  m.currarea,  m.talkstat,  m.compstat 

IF  SETCTALK")  =  "ON" 

SET  TALK  OFF 
m.UlksUt  =  "ON" 

ELSE 

m.talkstat  =  "OFF" 

ENDIF 

m.compstat  =  SET("COMPATIBLE") 

SET  COMPATIBLE  FOXPLUS 

m.rborder  =  SET("READBORDER") 

SET  READBORDER  ON 

m.cumrea  =  SELECTO 


•  ANALYSIS/Windows  Databases,  Indexes,  Relations 


IFUSED("class") 

SELECT  class 

SET  ORDER  TO  TAG  "class  id" 

ELSE 

SELECTO 

USE(LOCFILE("\smis\dbfs\class.dbr,"DBP,"Whcrc  is  class?")); 
AGAIN  ALIAS  class; 

ORDER  TAG  "class  id" 

ENDIF 


IF  USED("scction") 


SELECT  section 

SET  ORDER  TO  TAG  "sect  id" 

ELSE 

SELECT 0 

USE  (LOCFILE(’\5inis\db£s\5cction.dbr,*DBP, "Where  is  section?")); 
AGAIN  ALIAS  section ; 

ORDER  TAG  "sectjd" 

ENDIF 

IFUSED("eIemcnt") 

SELECT  element 

SET  ORDER  TO  TAG  "elm  id" 

ELSE 

SELECT  0 

USE  (LOCnLE("\smis\dbfs\element.dbr, "DBF, "Where  is  element?")); 
AGAIN  ALIAS  element ; 

ORDER  TAG  "eimjd" 

ENDIF 

IFUSED("steel') 

SELECT  steel 

SET  ORDER  TO  TAG  "steeljd" 

ELSE 

SELECT 0 

USE  (LOCnLE("\smis\dbfs\steel.dbr,"DBF", "Where  is  steel?")); 

AGAIN  ALIAS  steel ; 

ORDER  TAG  "steeljd" 

ENDIF 

IF  USED("life") 

SELECT  life 
SET  ORDER  TOO 

ELSE 

SELECT 0 

USE  (LOCnLE("\5mis\dbfs\life.dbr, "DBF, "Where  is  life?")); 

AGAIN  ALIAS  life  ; 

ORDER  0 

ENDIF 

IFUSED("momcurv") 

SELECT  momcurv 
SET  ORDER  TOO 

ELSE 

SELECT  0 

USE  (LOCFILE("\smis\dbfs\momcuiv.dbr, "DBF, "Where  is  momcurv?" 
AGAIN  ALIAS  momcurv ; 

ORDER  0 

ENDIF 


SELECT  class 


•  •  Windows  Window  deOnitions 

•  • 


IPNOTWEXISTCsmis'); 

OR  UPPER(WnTLE("SMIS*))  *»  'SMIS-PIX" ; 
OR  UPPER(WnTLE("SMlS"))  "SMIS.SCX*  ; 
OR  UPPER(WTITLE("SMIS*))  **  ’SMIS.MNX*  ; 
OR  UPPER(WnTLE("SMIS"))  ==  "SMIS-PRO" ; 
OR  UPPER(WTrrLE("SMIS")) ««  'SMIS.FRX" ; 
OR  UPPER(WnTLE(’SMIS*)) "SMIS.QPR" 
DEHNE  WINDOW  smis ; 

AT  3.083, 8.750; 

SIZE  24.538,100.200; 

TITLE  "Section  Analysis" ; 

FONT  "MS  Sans  Scrir,  8  ; 

FLOAT; 

NOCLOSE ; 

MINIMIZE ; 

SYSTEM 

ENDIF 


ANALYSIS/Windows  Setup  Code  -  SECTION  2 


#REGION  1 

PUBLIC  m.classjd,  m.sectjd,  m.dsJife,m.ntiines,m.ds_load,; 

m.loadbias,ni.loadcov 
DIMENSION  sections(40,2) 

SET  UDFPARMS  TO  REFERENCE 

STORE '  ’  TO  sections 

SELECT  DISTINCT  ALLTRIM(UPPER(cIass.classnanic)),  c]ass.class  id; 
FROM  CLASS; 

ORDER  BY  ciass.classname; 

INTO  ARRAY  classes 

in.betainin  s  0.0 
m.nlimes  =  0.0 

m.swbog  s  0.0 
in.b_swhog  =  1 .0 
m.covjswbog  a  0.2 
m.vwbog  a  0.0 
ni.b_vwbog  a  1 .0 
in.cov_vwbog  a  0.2 


m.b^swsag  *  1.0 
fn.cov_swsag  =  0.2 
m.vwsag  s  0.0 
m.b^vwsag  =1.0 
in.cov_vwsag  =  0.2 

in.session  = ' ' 


* 

•  ANALYSIS/Windows  Screen  Layout 


#REGIONl 

IFWVISIBLE("smis') 

ACTIVATE  WINDOW  smis  SAME 

ELSE 

ACTIVATE  WINDOW  smis  NOSHOW 

ENDIF 

@  3.769,57.600  SAY  "Minimum  BeU  level:"  +  CHR(13)  +  ; 
""  +  CHR(13)  +  ; 

"Time  Step  (in  years):"  ; 

SIZE  3.000,19.800,  0.000 ; 

FONT  "MS  Sans  SeriP,  8  ; 

STYLE  "T 

@  1.923,57.600  SAY  "SESSION  PARAMETERS"  ; 

FONT  "MS  Sans  Serif,  8  ; 

STYLE  "BP 

@  3.769,4.200  SAY  "Vessel  Class:"  ; 

FONT  "MS  Sans  Serif,  8  ; 

STYLE  "P 

@  5.615,4.200  SAY  "SecUon:*  ; 

FONT  "MS  Sans  Serif,  8  ; 

STYLE "P 

@  1.923,4.800  SAY  "SECTION  IDENTIHCATION"  ; 
FONT  "MS  Sans  Serif,  8  ; 

STYLE  "BP 

@  13.923,4.800  SAY  "Mean:*  ; 

FONT  *MS  Sans  Serif,  8  ; 

STYLE  "P 

@  15.769,4.800  SAY  "Bias:"  +  GHR(13)  +  ; 

*"  +  CHR(13)  +  ; 

•Coef.ofVar.:"  ; 

SIZE  3.000,13.000, 0.000 ; 

FONT  "MS  Sans  Serif,  8  ; 

STYLE  "P 

@9308,45.600  SAY  "Design  LOADS:"  ; 

FONT  "MS  Sans  Serif,  8  ; 

STYLE  "BP 


VWBM"  : 


@  8.308,2.600  TO  8.308,98.200  ; 

PEN  1,8; 

STYLET 

@  12.000,55.200  TO  19.462,55.200  ; 

PEN  1,8 

@  1.000,52.800  TO  8.385,52.800  ; 

PEN  1,  8 

@  19.385,3.000  TO  19.385,98.600  ; 
PEN  1,  8  ; 

STYLE  "1" 

@  12.615,23.600  SAY  "SWBM 

FONT  “MS  Sans  Serir,  8  ; 
STYLE  "T* 

@  12.615,64.800  SAY  “SWBM 

FONT  "MS  Sans  Serir,  8  ; 
STYLE T 

@  3.692,21.600  GET  m.which  class  ; 
PICTURE  ; 

FROM  classes  ; 

SIZE  1.538,22.167; 

DEFAULT  1 ; 

FONT  "MS  Sans  Serir,  8  ; 
STYLE "B"  ; 
WHEN_qldlcnsiO; 

VALID  _qld  IcfnivcO 
@  5.538,21.600  GET  in.wbicb_sect  ; 
PICTURE  ; 

FROM  sections  ; 

SIZE  1.538,22.167; 

DEFAULT  1 ; 

FONT  "MS  Sans  Serir,  8  ; 
STYLE "B"  ; 
WHEN_qldlcfni60; 

VALID  _qldlcfnnit0  ; 
DISABLE 

@  3.769,79.600  GET  m.betamin  ; 

SIZE  1.000,8.800 ; 

DEFAULT  0 ; 

FONT  "MS  Sans  Serir,  8  ; 
PICTURE  "@K  99.999" 

@  5.615,82.000  GET  m.timestep  ; 

SIZE  1.000,6.400 ; 

DEFAULT 0; 

FONT  "MS  Sans  Serir,  8  ; 
PICTURE  "@K" 

@  13.923,22.000  GET  m.swhog ; 

SIZE  1.000,13.600 ; 

DEFAULT 0; 

FONT  "MS  Sans  Serir,  8  ; 
PICTURE  "@KZ  999999999.9" 
DISABLE 

@  15.769,24.400  GET  in.b_swbog  ; 
SIZE  1.000,6.400” 

DEFAULT " " ; 


VWBM"  ; 


FONT  "MS  Sins  SeriP,  8  ; 
PICTURE ’@K' ; 

DISABLE 

@  17.615,24.400  GET  m.cov_swhog ; 
SIZE  1.000,6.400 :  ~ 
DEFAULT ■ • ; 

FONT  "MS  Sans  SeriP,  8  ; 
PICTURE  "@K" ; 

DISABLE 

@  13.92338.800  GET  m.vwhog  ; 

SIZE  1.000,13.600 : 
DEFAULT 0; 

FONT  "MS  Sins  Scrip,  8  ; 
PICTURE  "@KZ  99999999.9" 
DISABLE 

@  15.769,41.200  GET  in.b_vwhog  ; 
SIZE  1.000,6.400" 

DEFAULT  -  ■  ; 

FONT  "MS  Sans  SeriP,  8  ; 
PICTURE  "@K"  ; 

DISABLE 

@  17.615,41.200  GET  in.cov  vwhog  ; 
SIZE  1.000,6.400 ;  " 
DEFAULT ' ■ ; 

FONT  "MS  Sans  ScriP,  8  ; 
PICTURE  "@K"  ; 

DISABLE 

@  13.923,62.800  GET  m.swsag  ; 

SIZE  1.000,13.600 ; 
DEFAULT 0; 

FONT  'MS  Sans  SeriP,  8  ; 
PICTURE  "@KZ  999999999.9 
DISABLE 

@  15.769,65.200  GET  m.b  swsag  ; 
SIZE  1.000,6.400" 

DEFAULT "  • ; 

FONT  "MS  Sans  SeriP,  8  ; 
PICTURE  "@K"  ; 

DISABLE 

@  17.615,65300  GET  m.cov  swsag ; 
SIZE  1.000,6.400;  " 
DEFAULT""; 

FONT  "MS  Sans  SeriP,  8  ; 
PICTURE  "@K" ; 

DISABLE 

@  13.923,79.600  GET  m.vwsag ; 

SIZE  1.000,13.600 ; 
DEFAULT 0; 

FONT  "MS  Sans  SeriP,  8  ; 
PICTURE  "@KZ  999999999.9 
DISABLE 

@  15.769,82.000  GET  m.b_vwsag  ; 
SIZE  1.000,6.400" 

DEFAULT " " ; 


IFUSED("cliss*) 

SELECT  class 
USE 

ENDIF 

IF  USEDCscction") 

SELECT  section 
USE 

ENDIF 

IFUSEDCclemcnt’) 

SELECT  element 
USE 

ENDIF 

IF  USEDCsteel") 

SELECT  steel 
USE 

ENDIF 

IF  USED("life") 

SELECT  life 
USE 

ENDIF 

IF  USED("monKurv') 

SELECT  niomcurv 
USE 

ENDIF 

SELECT  (m.currarea) 


#REGION  0 

SET  READBORDER  &rborder 

IF  m.talkstat  =  ’ON” 

SET  TALK  ON 

ENDIF 

IF  m.compstat » "ON" 

SET  COMPATIBLE  ON 

ENDIF 


ANALYSIS/Windows  Cleanup  Code 


#REGION  1 


AND  Steel^teel_id  *  Element.steeijd; 
AND  Section^ct  id  >  nusect  id; 
INTO  ARRAY  temp 


m.count  s  ALEN(temp) 
n_elins  s  INT(m.count/16) 

m^ep  s  1 

DO  WHILE  m.bcU  >  m.BcUmin 

dimension  elms(n_eims,16) 

lastbeta  =  m.beta 

m.ycar  =  (m.step-l)*m.timestep 

DO  CORRODER 

DO  CAPACITY 

m.caphog  =  m.caphog/100CXXX) 
m.capsag  =  m.capsag/1000000 

m.beta  bog  =  (m.caphog  -  m.dmdbog); 

/sqrt((cov_capbog*m.capbog)''2  +  vdmdbog''2) 

m.betasag  =  (m.capsag  •  m.dmdsag); 

/sqrl((cov_capsag*m.capsag)''2  +  vdiudsag''2) 

m.beta  =  betabog 

SELECT  life 

APPEND  BLANK 
GATHER  MEMVAR 

m.step  =  m.slep  ■¥  1 


ENDDO 

m.step  «  m.step  - 1 

TTF  s  (m.step-l)*m.timestep  -f  (m.Betamin-lastbeta)  *; 
m.timestep/(m.beta-lastbeta) 

?TTF 

WATT  WINDOW  'YOU  WON' 

Set  up  wastage  limit  report 


elmsQ.lO)  =  temp(i-f  12) 
elmsQ.ll)  *  tcnip(i+13) 
elnis(j,12)  s  tcmp(i-f  14) 
eliie(j>^3)  “  teinp(i+15) 


i  »  i  +  16 
j.j  +  1 

IF  i  >  m.count 
done  «  .T. 
ENDIF 


ENDDO 
RETURN 
••••  END 


PROCEDURE  CAPACITY 


Input: 

m.count 

Output: 


elms  array  -  contains  the  element  specific  info  for  the  section 
•  length  of  elms  array 

Mmt_crv  -  Moment  vs  Curvature  array 


ailed  by:  ANALYSIS 
ails:  ELMPROPS 

STRSSTRN 


DIMENSION  el_curv(n_elms),  stiss(n_eims) 
Initialize... 


m.sect_area  »  0.0 
y_tiines_a  *  0.0 
I_nodes  s  0.0 
cr_stis  =  0.0 
I_node  =  0.0 
area  =  0.0 


Find  elastic  limit  curvature,  curvature  at  which  first  element  reaches 
its  elastic  limit,  stress  y_strs 

elaslim  s  10.  &&  impossibly  high 

y_stm  *  0.000000000000 
ef  curv  «  0.000000000000 

calculate  elastic  limit  curvature  in  the  hogging  condition 

FOR  i  =  1  TO  n^clras 

y_na  s  elms(i,9)  -  m.na 

y_stm  =  clms(i,12)/clms(i,10) 

IF  elms(i,14)  <  clms(i,12) 

IF  y_na  <  0 

el  curv  =  e!ms(i,14)/(elms(i,10)*y  na) 

ELSE 

el  curv  =  y  slm/y  na 

ENDIF 

ELSE 

el  curv  =  y  stm/y  na 

ENDIF 

elaslim  =  MIN(ABS(el  curv), elaslim) 

ENDFOR 

done  =  .F. 
j=2 

curvtr(j)  s  elaslim 
stress  =  0.0 

postj  =  5 


DO  WHILE  .NOT.  done 
k=j-i 

forcesum  =  0.0  &&  reset  values 

ea_sum  =  0.0 
niomcni(j)  =  0.0 

FOR  i  =  1  TO  n^clms 

strain  j=  curvtr(j)*(elms(i,9)  -  m.na) 

DO  stisstm  WITH  elms(i,10),elms(i,l  l),e!ms(i,12); 

,e]ms(i, 14), strain 

strss(i)  =  stress 

forcesum  =  forcesum  +  stress*clms(i,15) 
ea  sum  ^  ea_sum  ^  clms(i,10)^elms(i,15) 
ENDFOR  " 

IFj>2 

shift  -  forcesum/(ca_sum*curvtr(j)) 

ELSE 
shift  =  0.0 


y_na  *  elnis(i,9)  -  m.na 

y_strn  *  elms(i,12)/clins(i,10) 

IF  elms(i,14)  <  eiins(i,12) 

IF  y_M  >  0 

cl  curv  =  elms(i,14)/(elms(i,10)*y  na) 

ELSE 

el  curv  s  y  stm/y  na 

ENDIF 

ELSE 

el  curv  ®  y  stm/y  na 

ENDIF 

elaslim  =  -MIN(ABS(cl_curv),ABS(claslim)) 

ENDFOR 


done  s  .F. 
j=2 

cuMr(j)  =  elaslim 
stress  s  0.0 

forcesum  s  0.0 
ca^sum  =  0.0 
postj  s  5 

DO  WHILE  .NOT.  done 
k=j-i 

forcesum  =  0.0  &&  reset  values 

ea^sum  «  0.0 
moment(j) »  0.0 

FOR  i  s  1  TO  n_elms 

strain  =  curvtr(j)*(clms(i,9)  -  m.na) 

DO  strsstm  WITH  clms(i,10),elms(i,l  l),elms(i,12); 

»elms(i,14),strain 

strss(i)  s  stress 

forcesum  s  forcesum  stress*clms(i,15) 
ca  sum  sea  sum -f  elms(i,10)*elms(i,15) 
ENDFOR  " 


shift  s  forccsum/(ca_sum*curvtr(j)) 


m.na  s  m.na  shift 
moment(j)  s  0.0 
FOR  i  s  1  TO  n^clms 

moment(j)  s  moment(j)  +; 

strss(i)*elins(i»15)*(elms(i,9)*m.na) 


ENDFOR 


m.momnt  s  moment(j) 
m.curvatr  =  curvtr(j) 


in.u_strs,in.y_stn,m.poisson 
last  s  0.0 

Calculate  critical  stress  for  plate  between  stiffeners 

kD  ■  4*in.elasniod*m.pt''3/(12*(l.ni.poisson''2)) 
plt_cr  »  kD*PI0'‘2/(m.pfin.pb''2) 


done  « .F. 

kbeta  « (pb/pt)*sqft(u_strs/elasnKKl) 
pbe  «  m.pb 

area  «  in.pb*in.pt  m.wd*m.wt  -i-  in.fb*m.ft 

Calculate  column  (Euler)  buckling  stress 

e_area  «  m.pb*m.pt  +  m.wd*m.wt  +  ni.fb’m.ft 

centroid  *  (0.5*m.wd''2*in.wt  +  m.fb*ni.fl*m.wd)/e_area 

I_cent  a  (m.pb*m.pt''3  +  m.wt*m.wd''3  +  m.fb*m.ft''3)/12  +  ; 

m.pb*m.pt*cenlroid*2  +  m.wd*m.wt*(0.5*m.wd  -  centroid)'‘2; 
♦  m.fb*m.ft*(m.wd  •  centroid)''2 

gyradius  a  S0RT(I_cent/e_are8) 

col_cr  a  elasmod*(PI0*gyradius/m.space)''2 

Calculate  I_node 

crjstis  a  l,8*col_cr 

I_plate  a  (m.pt*m.pb''3/12)*SIN(m.theta)''2 
I_web  a  m.wt*m.wd''3*COS(m.tbeU)''2^ 

LOg  a  (m.ft*in.n)*3/12)*SIN(m.theU)*2 

l_node  a  I_plate  +  (I_web  +  m.wd*ni.wt*(Oii*m.wd*COS(m.theta))*2); 

+  (l_flg  +  m.fl*m.fb*(m.wd*COS(m.theta))''2) 


RETURN 


. . END 


Procedure  STRSSTRN 


CR_STRN  =  U_S'mN 

STRAIN  =  -STRAIN  &&.  one-sided  material  stress-strain  erv 

ELSE 

TENSION  =  .F. 

ENDIF 

IF  (STRAIN  >  CR^STRN) 

IF  (HARDSPOT) 

STRESS  =  m.U_STRS 

ELSE 

STRESS  »  CR_STRS  +  (CR  STRN  -  STRAIN)  •O.rm.ELASMOD 

ENDIF 

ELSE 

IF  (STRAIN  <*  Y_STRN) 

STRESS  =  STRAIN*in.ELASMOD . 

ELSE 

IF  (STRAIN  >  Y  STRN  jnd.  STRAIN  <  U_STRN) 

EXPONl  =  ^P((STRAIN  -  Y  STRN)^1) 

STRESS  =  (in.U_STRS/K2)*EXPONl/(l+EXPONl/k2) 

ELSE 

STRESS  =  m.U  STRS 
ENDIF 

ENDIF 

ENDIF 

IF  TENSION 

STRESS  » -STRESS 


ENDIF 

* 

RETURN 

*  « 

•  •  _QLD1CFLSI 

«  * 

ni.which_class  WHEN 

*  *  Function  Origin: 

•  •  From  Platfomt: 

Windows 

•  •  From  Screen: 

ANALYSIS,  Record  Number:  22 

•  •  Variable: 

m.whicb^class 

•  •  Called  By: 

WHEN  Clause 

•  •  Snippet  Number: 

•  • 

1 

FUNCTION  qldlcHsi 

&&  m.which_class  WHEN 

#REGION  1  ” 

m.wlucb_class  s  1 

_QLD1CFMVC  m.which_class  VALID 
Function  Origin: 


WHERE  scctionxlassjd  =  m.classjd; 
ORDER  BY  sectjd; 

INTO  ARRAY  sections 

SHOW  GET  in.which.secl  ENABLE 
.CUROBJ  =  OBJNUM(which.sccl) 

*  find  the  correct  class  record  and  point  to  it 


.QLD1CFNI6  in.which.sect  WHEN 
Function  Origin: 

From  Plalfomi:  Windows 

From  Screen:  ANALYSIS,  Record  Number:  2S 

Variable:  m.which  sect 

Called  By:  WHEN  Clause 

Snippet  Number:  3 


FUNCTION  qldlcfnib  SlSl  m.which  sect  WHEN 
#REGION  1  " 
m.which.secl  =  1 


.QLDICFNMT  m.which.secl  VALID 

Function  Origin: 

From  Platfonn:  Windows 

From  Screen:  ANALYSIS,  Record  Number:  23 

Variable:  m.which.sect 

Called  By:  VALIDQause 

Snippet  Number:  4 


FUNCTION  qldlcfnmt  &&  m.which  sect  VALID 
#REGION  1  “ 

m.sectjd  =  sections(m.which_sect,2) 


.QLDICFOOM  m.begin  VALID 
Function  Origin: 


Appendix  B:  Section  Ideauzation  Table 


ELM  ID 

PB  PT 

PR 

WD 

WT 

WR 

FB 

FT 

FR 

Y  NODE 

THETA 

STEEL  ID 

BOIP 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.126 

0.00 

2222 

B02P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.126 

0.00 

2222 

B03P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B04P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B05P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BOSS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B06P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BOSS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B07P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B07S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B08P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BOBS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B09P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B09S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BIOP 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BIOS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

Blip 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BUS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B12P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B12S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B13P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B13S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B14P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B14S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B15P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BISS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B16P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B16S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B17P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B17S 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B18P 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

BIBS 

1000.0  2S.0 

0.197 

797.0 

16.0 

0.063 

200.0 

33.0 

0.063 

0.125 

0.00 

2222 

B19P 

9SO.O  18.0 

0.061 

297.0 

11.6 

0.036 

100.0 

16.0 

0.060 

950.000 

1.07 

2222 

B19S 

9S0.0  18.0 

0.061 

297.0 

11.6 

0.036 

100.0 

16.0 

0.060 

950.000 

1.07 

2222 

B20P 

42S.0  2S.0 

0.036 

476.0 

18.0 

0.061 

0.0 

0.0 

0.000 

1900.000 

3.14 

2222 

B20S 

42S.0  2S.0 

0.036 

476.0 

18.0 

0.061 

0.0 

0.0 

0.000 

1900.000 

3.14 

2222 

B21P 

9S0.0  18.0 

0.061 

297.0 

11.6 

0.036 

100.0 

16.0 

0.060 

950.000 

1.07 

2222 

B21S 

9S0.0  18.0 

0.061 

297.0 

11.6 

0.036 

100.0 

16.0 

0.060 

950.000 

1.07 

2222 

B22P 

370.0  16.0 

0.036 

476.0 

18.0 

0.061 

0.0 

0.0 

0.000 

1900.000 

3.14 

2222 

B22S 

370.0  16.0 

0.036 

476.0 

18.0 

0.061 

0.0 

0.0 

0.000 

1900.000 

3.14 

2222 

D01P 

1000.0  2S.0 

0.110 

480.0 

32.0 

0.063 

0.0 

0.0 

0.000 

25800.000 

3.14 

2222 

DOTS 

1000.0  26.0 

0.110 

480.0 

32.0 

0.063 

0.0 

0.0 

0.000 

25800.000 

3.14 

2222 

D02P 

1000.0  26.0 

0.110 

480.0 

32.0 

0.063 

0.0 

0.0 

0.000 

25800.000 

3.14 

2222 
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